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SECTION 1 


INTRODUCTION 


1. 1 BACKGROUND 

Work under two concurrent NASA contracts (NASS -21453 and NAS8-30503) 
resulted in publication of the NASA Handbook for Estimating Toxic Fuel Hazards 
(Dumbauld, Bjorklund, Cramer and Record, 1970). The handbook was accompanied 
by a computer program (NASA/MSFC Multilayer Diffusion Models Program) specifi- 
cally designed for research -planning estimates of downwind dispersions from normal 
and abnormal launches of rocket vehicles as well as from accidental cold spills and 
leaks of toxic fuels. The transport and diffusion models that formed the basis of 
the program were extensions of the generalized concentration, dosage and deposition 
model concept originally developed for the U. S. Army (Cramer, et al. , 1964; 1967) 
to include a multilayer tropospheric construct with provision for changes in atmospheric 
structure along the cloud trajectory downwind from the launch site. As experience 
was gained in the application of the NASA/MSFC Multilayer Diffusion Models Program, 
revisions were made to simplify its use and to incorporate the latest available 
advances in diffusion technology. Version 2 of the Program, designed for use in 
operational support of vehicle launches, was completed in 1973. A technical report 
by Dumbauld, Bjorklund and Bowers (1973) describes the updated transport, diffusion 
and cloud rise models and the operation of the computer program. 

Since 1973, work has continued on the development and refinement of the 
models and the computer program. 

1. 2 PURPOSE 

The purpose of this report is to document; 

• The latest refinements and developments in the transport and 
diffusion models 
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• Version 5 of the NASA/MSFC Multilayer Diffusion Models Program 

• The Preprocessor Program used in conjunction with Version 5 of the 
Multilayer Program 

The majorify of the changes in the NASA/MSFC Multilayer Diffusion Models Program 
have been made to further simplify the use of the models and computer program so 
that users can apply the techniques for hazard estimation with only a minimum 
knowledge of the technical aspects of diffusion meteorology and associated disciplines. 
The simplification of program use was accomplished in two ways. First, a Prepro- 
cessor Program was designed to automate the calculation of the: 

• Effective height of the stabilized exhaust cloud 

• Distribution of exhaust products in the meteorological layers 
containing the cloud 

• Meteorological inputs to the transport and diffusion models 

Use of the Preprocessor Program essentially requires only surface meteorological 
and radiosonde data be available for input. Second, the NASA/MSFC Multilayer 
Di^sion Models Program was redesigned for more efficient operation and an 
extensive capability for automated plotting of isopleths and centerline profiles of 
concentration, dosage and deposition due to gravitational settling and precipitation 
scavenging is now provided to assist users in interpretation of computational results. 

In addition to simplifications in the use of computer programs, technical 
features of the transport and diffusion models have been altered to more accurately 
reflect observations of cloud rise during actual launches of solid-fueled vehicles 
and the physical aspects of the transport and diffusion process. For example, altera- 
tions have been made in assumptions regarding the physical shape of the ground 
cloud at time of cloud stabilization and in the diffusion models to account for partial 
reflection (or partial absorption) of material at air-surface interfaces. 
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1.3 


ORGANIZATION OF THE REPORT 


The main body of the report contaias five sectiODS. Section 2 contains a 
description of the models and algorithms contained in the Preprocessor Program. 

The generalized transport and diffusion models for calculating concentration, 
dosage and deposition contained in the NASA/MSFC Multilayer Diffusion Models 
Program are described in Section 3, A brief description of the Preprocessor and 
Diffusion Models Programs is given in Section 4 and Section 5 contains some 
sample problems and their solutions obtained using the Programs. 

There are five appendices to the report. Appendix A contains user instruc- 
tions for implementing the Preprocessor Program and Appendix B contains user 
instructions for implementing the NASA/MSFC Multilayer Diffusion Models Program. 
Computer program listings for both the progprams are given in Appendix C. Appendix 
D contains example computer output listings for the sample problems described In 
Section 5 of the main text. Finally, a derivation of the vertical term for Model 4 
of the NASA/MSFC Multilayer Diffusion Models Program is presented in Appendix 
E. 
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SECTION 2 


PREPROCESSOR MODELS AND ALGORITHMS 

A Preprocessor Program has beea developed to simplify the calculation of 
source and meteorological inputs required by the NASA/MSFC Multilayer Diffusion 
Models Program. The models and algorithms incorporated in the Preprocessor 
Program for calculating the requisite model inputs are described in this section of 
the report. In its present form, the Preprocessor Program has the capability of 
processing data for the Space Shuttle, Titan III, Delta-Thor and Minuteman II 
vehicles. 


2. 1 FUEL PROPERTIES AND VEHICLE RISE DATA 

Properties of the vehicle fuel and rise data are required to calculate cloud 
rise and the source strength distribution of atmospheric pollutants in the troposphere 
resulting from normal and abnormal launches. Two types of abnormal launches 
have been hypothesized for the Space Shuttle and Titan III vehicles. In one situation, 
it is assumed that a single solid engine of the Titan or Space Shuttle zero stage 
ignites and bums over the normal engine firing period while the vehicle remains in 
a hold-down configuration on the pad (single-engine burn). In the second pad -abort 
situation, it is assumed that an on-pad explosion ruptures the casings of the two 
solid engines and that all the solid propellant then falls to the ground in the vicinity 
of the pad and bums at a constant rate over a 5 -minute period (slow burn). In both 
pad-abort situations, it is further assumed that the other vehicle stages are unaffected 
by the burning of the zero -stage solid propellant and do not therefore contribute in 
any way to the combustion products or heat releasedducing-.the onrpad, aborts. .. The 
net affects of this latter assumption is to minimize the cloud rise and maximize the 
concentration of the pollutants produced by the burning of the solid propellant. Slow- 
burn incidents in which fuel burns over a 5-minute period are also hypothesized for 
the Delta-Thor and Minuteman II vehicles. In the case of the Delta-Thor, both the 
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solid propellant from the six castor engines and liquid fuel from the zero-stage 
engine are assumed to contribute to the heat available for cloud rise and to pollutant 
concentrations. The solid propellant from all three stages of the Minuteman II is 
assumed to fall to the ground near the launch site as the vehicle is destroyed just after 
clearing the launch silo. 

Fuel expenditure and heat content data incorporated in the Preprocessor 
Program for normal and abnormal launches of the four types of vehicles are given 
in Table 2-1. Fuel expenditure rates for normal launches were obtained by averaging 
fuel consumption over the approximate period from lift-off until the vehicle is about 
3 kilometers above the surface. For single -engine burns, the fuel expenditure rate 
was calculated by dividing the total propellant weight in a single engine by the normal 
firing period shown in the table. Similarly for slow burns in the vicinity of the 
launch area', the fuel expenditure rate was obtained by dividing the total weight of 
all solid engines, and the zero stage liquid engine in case of the Delta-Thor, by the 
total assumed bum time of 300 seconds. 

The fuel heat contents used in calculating cloud rise arc estimates based 
on the best available information and are subject to change as additional information 
becomes available. In case of the Tital III vehicle, the fuel heat content is based 
on estimates from two-phase (gas and solid) flow accounting approximately for heat 
gains from afterburning and heat losses due to radiation, and on the experience to 
date in predicting cloud rise for a number of actual launches using the instantaneous 
cloud rise model described in Section 2. 2 below. Lacking better information, the 
same heat content has also been used for the solid engines of the Space Shuttle vehicle. 
A heat content value of 500 calories per gram was used for the liquid fuel in the zero 
stage engines of the Space Shuttle and Delta-Thor vehicles because previous experience 
in predicting cloud rise from Saturn launches indicates this value may be appropriate. 
The heat contents of the solid fuel for the Delta-Thor and Minuteman II are based on 
early estimates of the heat content for solid fueled engines and are probably low. 
However, experience in predicting cloud rise from a limited number of Delta-Thor 
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TABLE 2-1 


PUEL EXPENDITURE AND HEAT CONTENT DATA 



Vehicle Type 

Property 

Space Shuttle 

Titan III 

Delta -Thor 

Minuteman It 


(a) 

Normal Launch 


Fuel Expenditure 
Rate (g sec”^) 





Solid Engine 

9. 335 X 106 

4. 174 X 10® 

6. 05 X 105 

3.771 xlO^ 

Liquid Engine 

1. 531 X 10® 


3. 13 X 10® 


Effective Fuel Heat 





Content (cal g~^) 





Solid Fuel 

2500 

2500 

691 

691 

Liquid Fuel 

500 


500 



(b) Single Engine Burn 


Fuel Expenditure 
Rate (g sec~^) 

3.753 x 106 

1. 742 X 10® 



Normal Firing 
Period (sec) 

122 

122 



Effective Fuel Heat 





Content (cal g~^) 

1274 

1036 
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TABLE 2-1 


FUEL EXPENDITURE AND HEAT CONTENT DATA 
(CONTINUED) 


Property 


Vehicle 

1 

} Type 

Space Shuttle 

Titaa III 

Delta-Thor 

Mlauteman n 


(c) Slow Bum 


Fuel Expenditure 
Rate (g sec~^) 





Solid Engine 
Liquid Engine 

3.052 X 10® 

1. 301 X 10® 

7.462 X 10^ 
2.212 X 10® 

2. 87 X 10*^ 

Total Burn Time 
(sec) 

300 

300 

300 

300 

Effective Fuel Heat 
Content (cal g"^) 

1000 

1000 

1000 

1000 















launches using the heat contents in Table 2-1 and a combination of the instantaneous 
and continuous cloud rise models indicate serious errors are not being made in the 
rise prediction. In any case, the conservative low estimates of heat available for cloud 
rise for the Delta-Thor and Minuteman II vehicles have the effect of minimizing 
cloud rise and maximizing ground-level concentrations. There is no experience in 
predicting cloud rise from launch aborts of any of the four vehicle types. The heat 
contents shown in Table 2-1 for single-engine burns of the Space Shuttle and Titan 
vehicles are based on the estimate that about 1500 calories per gram are available 
from the fuel burn from a single engine, but that some of the available heat is dissi- 

3 

pated in heating and vaporizing 1. 26 x 10 kilograms per second of water used to 
spray the launch pad during the burn. The heat content of 1000 calories per gram 
hypothesized for slow burns of fuel from the zero-stages of the vehicle is thought to . 
be a realistic estimate of the heat available for cloud rise from the burning of uncon- 
fined fuel. 

The fraction by weight of pollutants comprising the rocket exhaust products of 
the four tjrpes of vehicles used in the calculation of the vertical distribution of pollu- 
tants in the lower troposphere are shown in Table 2-2. Two sets of fractions are 
given for the Minuteman II because in an abnormal launch all three stages are hypoth- 
esized as being destroyed and in normal launches only the first stage contributes to 
the pollutant distribution in the lower troposphere. Recent indications (Cicerone, 
Stedman and Stolarski, 1973) are that CO in the exhaust plume may quickly oxidize 
to C02» in which case no toxic CO problem exists. However, because the oxidation 
of CO to CO 2 has not been verified by direct measurements, the distribution of CO 
in the lower troposphere is calculated in the Preprocessor Program. 

It should be noted that the values of heat content and fuel expenditure rates in 
Tabel 2-1 and the fraction of pollutants by weight in Table 2-2 may differ from the 
values assigned as constants in the Preprocessor Program described in Appendix A. 

This occurs because, in the Rrogram, the constants are manipulated in different ways 
to calculate cloud rise and the vertical distribution of pollutants. Thus, after the 
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TABLE 2-2 


POLLUTANT COMPOSITION OF FUEL (Fraction by Weight) 



Vehicle Type 

Pollutant 

Space Shuttle 

Titan III 

Dclta-Thor 

Minuteman n 





Normal 

Abnormal 

HCl 

0.207 

0.210 

0.208 

0.197 

0.204 

* 

CO 

Solid Engine 
Liquid Engine 

0.280 

0.279 

0.223 

0.473 

0.220 

0.219 

CO2 


0.029 




AI2O3 

0. 304 

0.304 

0.378 

0.277 

0.280 


May be oxidized to CO2 


TABLE 2-3 
VEHICLE RISE DATA 


Power-Law 

Vehicle Type 

Coefficients 

Space Shuttle 

Titan III 

Delta -Thor 

Minuteman II 

a 

0.664 

0.635 

1.321 

0.440 

b 

0.485 

0.484 

0.395 

0.479 
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calculations are completed by the Preprocessor Program, the values shown in Tables 
2-1 and 2-2 have been, in effect, used in the calculations. 

The altitude-time curves of the various vehicles are also required to calculate 
the buoyant rise of the ground cloud of exhaust products. A simple power-law rela- 
tionship is sufficiently accurate to describe the altitude -time curve in the first several 
thousand meters near the surface. A logarithmic least-squares regression curve of 
the form 

Time = a (Altitude)^ (2-1) 

where time is in seconds and altitude in meters was fitted to the altitude-time infor- 
mation for all vehicles. The resulting values of the coefficients a and b obtained from 
the fitting procedure are given in Table 2-3. 

2.2 CLOUD RISE FORMULAS 

The burning of rocket engines during launches and on-pad aborts results in 
the formation of a cloud of hot exhaust products which subsequently rises and entrains 
ambient air until an equilibrium with ambient conditions is reached. For normal 
launches, the cloud is formed principally by the forced ascent of hot, turbulent 
exhaust products that have been deflected laterally and vertically by the launch pad 
hardware and the ground surface. In the case of normal launches of solid-fueled 
vehicles or vehicles with a number of solid boosters, vehicle hold-down times are 
minimal and vehicle residence times in the lowest kilometer of the atmosphere are 
relatively short. The exhaust products contained in the stabilized ground cloud are 
therefore emitted over a time period of from about 10 to 30 seconds. Experience to 
date shows that the buoyant rise from vehicles with solid-fueled first stages is best 
predicted by using a cloud-rise model for instantaneous sources and rise from 
vehicles having liquid-fueled first stages is best predicted using a cloud-rise model 
for continuous sources. Thus, an instantaneous source cloud -rise model is used 
in the Preprocessor Program to predict buoyant rise of the ground cloud for normal 
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launches of the Space Shuttle, Titan III and Minuteman II vehicles. Limited experi- 
ence in predicting the buoyant rise of the ground cloud generated by normal launches 
of the Delta-Thor vehicle with its large liquid-fueled first stage and six solid-fueled 
boosters indicates that an average of the rise predicted using the instantaneous and 
continuous source cloud rise models is appropriate. While no cloud rise data are 
available for on-pad aborts of any of the four vehicle types, cloud rise data from 
static-tests of liquid fueled rocket engines indicate that the use of a cloud rise model 
for continuous sources is appropriate in these cases. 


The instantaneous and continuous cloud-rise models used in the Preprocessor 
Program are based on work by Briggs (1969, 1970). Derivations of these models 
also appear in the report by Dumbauld, _et al. (1973). Only the simplified forms of 
the models used in the Preprocessor Program are given below. 

2. 2. 1 Instantaneous Source Cloud-Rise Model 


The maximum cloud rise z jj^ from an instantaneous source in a thermally 
stable atmosphere is given by the expression 


z 


ml 


where 


F 


I 




initial buoyancy term 

3 g Q I 
4 Cp T 7T p 

gravitational acceleration (9.8 m sec ) 
effective heat released (cal) 
specific heat of air at constant pressure 
(cal g“^ ®K) 

ambient air temperature (°K) at the surface 
reference height Zjj 
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p 

s 


”3 

density of ambient air (g m ) at Zg 

S 9$ ^ g A* 

T 9z T Az 


A^ 

•Az 


vertical gradient of ambient virtual potential 
temperature 


n 


= entrainment constant - 0. 64 


The time t^j^ required for the center of the ground cloud to reach an altitude 
given by the expression 


t„ = s arccos 
Ik 




where 


pc A 4» 

P ^ A z 

* 

The time t required for the cloud to reach the stabilization height Zj^| is thus 
given by the egression 


t 


* 


Tt 8 


- 1/2 


In calculating z from Equation (2-2) the effective heat released for 
cloud rise Q| is calculated from the relationship 


Ql « Qp tR •j^ml} 

where 

Qj* = rate of neat released by the burning fuel 
= H • W 

H = fuel heat content 


is 

(2-3) 


(2-4) 

(2-5) 
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w 


fuel expenditure rate 


seconds for the vehicle to 
reach the height of the centroid of the 
stabilized ground cloud Zj^^j 

= ( 2 - 6 ) 

where a and b are the power-law coefficients from Table 2-3. 


According to the above formulas, the following quantities are interrelated: 
the calculated maximum cloud rise Zj^j, the height over which the virtual potential 
temperature gradient A$ / Az is measured, and the value of used in ob- 

taining Qj. The final value of maximum cloud rise must therefore be found through 
iteration of Equation (2-1). The Preprocessor Program accepts radiosonde data 
which provides the temperature, pressure and relative humidity profiles as ftmctions 
of height in the troposphere and the surface density p for use in the cloud rise 
calculation. The virtual potential temperature at each radiosonde measurement 
height Zj^ is obtained from the expression (Tabata, 1973) 



4*, = T, 

1+1.61 wj^ 


1000 

where 

k k 

1 + w, 

L k _ 


_^k _ 


T 


k 


Air temperature at z in degrees absolute 


Pj^ = barometric pressure at z (mb) 

= mixing ratio at zj^ 

0. 622 (RH) ® s ; k 
k ®s ; k 

(Rll)j^ = relative humidity at Zj^ expressed as a fraction 
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®s; k 

X = 

and c, d and e are constants given by 

c = 

d 

e = 


saturation vapor pressure at Zj^ 
^^{c-dx-ex^) 

1000/T, 

k 

8.42926604 
1. 82717843 
0. 071208271 


The iteration process then begins by assuming that the first level above the surface 
reference height (z^^ (k=l} at which a radiosonde measurement is avail- 
able (Zj^ {k = 2 } ) is equal to and solving Equation (2-6) for t^ j a 

and Equation (2-5) for Q^. The lapse rate of virtual potential temperature between 
the surface reference height and the first height above the surface is then obtained 
from the e^qjression 


A ® ^ {k = 2} - 4-jJk = 1} 

A z Zj^ {k = 2} - Zj^ 

and Equation (2-2) solved for z If the value of z^^^j calculated from Equation 
(2-6) is less than the value z^ {k = 2} , then 10 meters is subtracted from z ^ {k = 2} 
and the process repeated using the same value of A4/Az until the estimated value of 
Zmi is within i 10 meters of the value of Zj^^j calculated from Equation (2-2). If on 
the first iteration step the value of Zj^^j calculated from Equation (2-6) exceeds the 
value Zj^{k = 2}, then the value of zj^ {k = 3} is used as an estimate of z^^j and the 
iteration process is continued by increasing k until the value of Zj^^j from Equation 
(2-6) exceeds the estimated value of Zjj^j. However, for each iteration where k > 2, 
a least -square regression fit of the form 
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(2-7) 



is used, at the suggestion of Dr. Briscoe Stephens of NASA/MSFC, to obtain an 
estimate of the vertical potential temperature gradient. As a result of this iteration 
procedure, the Preprocessor Program obtains an estimated final height of the cloud 
centroid within 1 10 meters of the exact value that should be calculated. 

It should be noted that Equation (2-2) is for use when the atmosphere is 
thermally stable. For this reason, if the virtual potential temperature gradient is 
ever less than 3. 322 x 10“^ degrees per meter, the gradient is set equal to 3. 322 x 
10”^ degrees per meter. Experience has shown that use of this device yields nearly 
the same final cloud rise as the cloud rise that would have been obtained had a formula 
for an adiabatic atmosphere been used in the calculation. 


2. 2. 2 Continuous Source Cloud Bise Model 

The maximum cloud rise from a continuous source is 


where 


z 

me 



Fg = buoyancy flux 

. g 

.pOpT 

Qg = effective heat rate (cal sec ^) 


(2-8) 
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H* W 


u 

c 


u 


k 


Ui 


K 


eatrainmeat constant = 0. 5 

height-weighted mean wind speed 
between the surface and z 

me 

K 



index of the last radiosonde 
measurement point below the 
height of the cloud centroid at 
stabilization z 

me 


and Uj^ are wind speeds at the radiosonde observation heights Zy_. The time ^5k 
required for the centroid of the ground cloud to reach an altitude for continuous 
sources is given by the expression 



( 2 - 9 ) 


where 


J = 


^k 


p c Try 

p 


A ^ 
A z 


i = 1 
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The time t required for the cloud to reach the stabilization height z is identical 
to the time for instantaneous sources given by Equation (2-4). 

As in the case of instantaneous sources, iteration is required to solve for 
the height because it is interrelated with the height over which A VAz is 
measured and with the wind speed u^. The same iteration procedure described in 
Section 2.2. 1 above for instantaneous sources is used in the Preprocessor routine 
for calculating 

2. 3 CALCULATION OF SOURCE MODEL INPUT PARAMETERS 

The Preprocessor Program is used to calculate the dimensions and spatial 
position of the stabilized ground cloud and the distribution of exhaust products within 
the cloud. In its present form, the Preprocessor Program can be used only with 
Models 3 and 4 (see Section 3) to calculate the source model input parameters for 
those portions of the stabilized ground cloud that are- contained within the surface 
mixing layer. The source geometry identified \vith Models 3 and 4 is described in 
Section 5. In brief. Model 3 assumes that the portion of the stabilized ground cloud 
in the surface mixing layer is all contained within a spherical volume source. In 
Model 4, the stabilized ground cloud in the surface mixing layer is assumed to be 
contained in a number of cylindrical volume sources extending from the ground 
surface to the top of the mixing layer. Extension of the Preprocessor Program to 
include calculations of source inputs for all the models in the NASA/MSFC Multi- 
layer Diffusion Models Program is currently in process with emphasis on implemen- 
tation for the RE EDA system of NASA/MSFC. 

2. 3. 1 Calculation of the Dimensions and Spatial Position 

of the Stabilized Ground Cloud 

In the Preprocessor Program, the cloud radius at any 
height z during cloud rise is calculated by the expression 
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r{z} 


( 2 - 10 ) 


yz 


z < z 


m 


y (2z -z); z <z< 2 z 

m m ~ mi 


z > 2 z 


m 


where y is equal 7 ^ or y and z^^^ is equal to z^^j or z^^^ depending on whether the 
source is instantaneous (I) or continuous (c). 


The alongwind, crosswind and vertical source dimensions in the surface 
mixing layer for Model 3 are calculated under the following assumptions: 

• The alongwind, crosswind and vertical distributions of exhaust products 
within the stabilized cloud are Gaussian 

• The concentration of exhaust products at a distance of one radius from 
the cloud centroid is 10 percent of the concentration at the centroid 

Thus the standard deviations of the alongwind (o (K = 1}) , crosswind (o (K= 1}) 

xo y® 

and vertical (cr^^(K=l}) distribution which define the source dimensions at the point 
of cloud stabilization are calculated from the relationships 


a (K=l} = o- {K = l} = r{z } / 2. 15 = yz / 2. 15 
xo ‘ ' yo ' ' m' m 


( 2 - 11 ) 


o 

zo 


(K=l}= 


r{z }/2.15 
' m' 

z„^ (K=l}-z 
TK ' m 

4,3 




r (z } 

z -r(z )<z_„(K=l)<z + 

m ' m' TK' m 



( 2 - 12 ) 
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where 

z {K = 1} = depth of the surface mixing layer for Model 3 

TIv 


In the special case where the bulk of the cloud is above the surface mixing layer 

< the Preprocessor Program prints a message indicating 

that calculations are terminated for Model 3 since its use is invalid in the surface 
mixing layer for this case. The effective height {K = 1} of the stabilized cloud 
centroid in the surface mixing layer for Model 3 is given by 


m 


;z_„{K=l}>z +r{z } 
’ TK ' - m ' m' 


Hjj{K=l} = 


.(2-13) 


z_..{K=l}+z -r{z } 
— — ^ ^ 


In the case of Model 4, the Preprocessor Program assumes that the surface 
mixing layer, for which the depth is defined as = l) » is comprised of one or 

more sublayers K with boundaries at heights in the surface mixing layer specified 
in the radiosonde message. Alongwind and crosswind source dimensions are defined 
in each K sublayer for normal launches by 


yz 
' m 

2. 15 


» z’ < z 
— m 


<r 

xo 


(K}=o 


yo 


{K} 


y(2z -z') 
m 

2.15 

200 = 93 

2. 15 


V. 

and for abnormal launches by 


; z' >z , yz’ > 200 
m ' 

; z' >z , yz' < 200 
m ’ — 


(2-14) 
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m 

2.15 


» z' < z 
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cr {K} = cr {K} = -{ 
xo' ' yo ' ' 


y (2z - z') 

; z < z' < z 

2.15 m - m 


; z' > 2z 


m 


th 


midpoint of the K sublayer 




"TK 


th 

height of the top of the K sublayer 


BK 


th 

height of the base of the K sublayer 


The corresponding vertical source dimension for Model 4 is zero since the model 
assumes .a rectangular source distribution in the sublayers and uses error functions 
to simulate a vertical line source between sublayer boundaries which is comprised 
of an infinite number of (vertical) point sources (see Section 3.4 below). 

For Model 4, the spatial position in the horizontal plane of the cloud in any 
sublayer K at t*, the time of cloud stabilization, with respect to the origin at the 
laimch pad is assumed to be given by 




90 - tan”^ 

2 2 ^^ 


where 
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(2-15) 

(2-16) 

(2-17) 

(2-18) 
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For 


(2-19) 


t. , , or t , , < t 

I, k + 1 c, k + 1 






(t 1 ■ I,) 1 ■*■ ^*1,) / ^ » instantaneous source 
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t 1 “ t , ) (z. , - z.) (u. , + u.) / 2 

!, k + 1 c, k' ^ ' i + 1 i' ' 1 + 1 r 


i = 1 


continuous 
* source 


S ^^i + i ■ ^i^ 
i= 1 


( 2 - 20 ) 


where Q . and Uj^ are respectively the wind direction and wind speed from radiosonde 

measurements at the k height and tx i and t , are the times respectively from 

* th 

Equations (2-3) and (2-9) for the cloud to pass through the k height. When 

t- , , or t , , > t* ; 

I, k + 1 < 5 , k + 1 


(0 + e, ) / 2 

m k 


(2-21) 
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(2-23) 
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and 


R 

k 


.where 



R 

c 


R.. (t - t, , , ) ; instantaneous source 
I ' I,k+r 


f * 

R continuous source 

c c, k+1 


(2-24) 
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'z - z, 

m k 


( 


(2-26) 


i = 1 


For Model 3, only the last position of the stabilized cloud, obtained by stepping 
throu^ Equations (2-15) to (2-26), is used to identify the position of the cloud at 
stabilization* 


2. 3. 2 Calculation of the Distribution of Exhaust Products 

within the Stabilized Cloud 


The fraction by weight of pollutant material in the stabilized 
ground cloud in each of the K sublayers within the surface mixing layer for Model 4 
is given by the expression 


F (k} = 




(2-27) 


where 

Q 


total weight of exhaust products in the stabilized 
ground cloud 
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(2-28) 

FM 

= 

fraction by weight of pollutant material in the fuel 



from Table 2-2 
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Inspection of Equations (2 -29) and (2-30) shows that a Gaussian vertical distiibution 
of material is assumed about the height of the stabilized ground cloud z^. If the 
height z is less than the depth of the surface mixing layer z {L = 1} and the launch 
is normal, the vehicle exhaust trail will insert more material into the surface mixing 
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layer than given by Equation (2-27). In this case, the total fraction of material 
th 

in the K sublayer is given by 


FjKHW. <z^ (L=l> 


(K) . 


(2-31) 




In case of Model 3, only the total source strength of the pollutant in the surface 
mixing layer {K = 1} is required. Thus, 


F^iK=l} = QP{z^j^{K = l}} 


(2-32) 


where the top of the surface mixing layer, Q is defined by Equation (2-28) 

above and P{z^^(K=l}} is defined by Equation (2-29). 

Since the desired concentration units for HCl, CO and CO 2 are parts per 
million, the complete expression for Q„ is 


Q 


K 




^10^ mg^ 


S 


\f 22 . 4 \f T \/l013.2\ 
i )(-. (2-33) 

/ \ M /\273.16/ \ P / 


where 

M = molecular weight of HCl, CO, or CO^ 

T = ambient air temperature (°K) 

ambient pressure (mb) 


P 



For AI 2 O 3 , the desired concentration units are milligrams per cubic meter, so that 


Q 


K 


F^(K}(10^mg/g) 


(2-34) 


2. 4 CALCULATION OF METEOROLOGICAL MODEL INPUT PARAMETERS 


The majority of meteorological model input parameters required by the 
Preprocessor Rrogram, and by Models 3 and 4 described in Section 3 below are obtained 
from radiosonde observations. These include: 



wind speed (m sec~^) 

wind direction (deg) 
ambient air temperature (°K) 

ambient relative humidity expressed 
as a percentage 


Pj^ = atmospheric pressure (mb) 

“3 

p = surface air density (gm ) 


In addition, the use of Models 3 and 4 to predict concentrations and dosages in the 
suirface mixing layer requires the following meteorological parameters: 
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{L=l) or cr ^ {K=l} = mean layer standard deviation of the wind 

azimuth angle in radians for the surface 
mixing layer for Model 4 (L notation) and 
Model 3 (K notation). 

ff' t{L= 1} or cr' „{K=l} = mean layer standard deviation of the wind 

elevation angle in radians for the surface 
mixing layer for Model 4 and Model 3. 

a or Q! =t lateral diffusion coefficient in the surface 

mixing layer for Model 4 and Model 3. 


/3 _ or |S = vertical diffusion coefficient in the surface 

Lf K 

mixing layer for Model 4 and Mcxlel 3. 

z^j^{L=l) or z^j^ {K=l} = depth of the surface mixing layer for 

. Model 4 and Model 3. 


The standard deviation of the wind azimuth angle in the surface mixing 
layer, in the absence of turbulence induced in the layer by the rocket vehicle itself, 
is proportional to the source function time t and the wind speed profile in the mixing 
layer. Dumbauld, et al. (1970) give the following relationships: 


and 



T < T < 600 seconds 
o 




z 


z 


R 


•P 


(2-35) 


(2-36) 
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where 



P 


standard deviation of the wind aximuth 

angle for the reference time t measured 

o 

at a reference height z 

R 

power-law coefficient of the wind speed 
profile in the surface mixing layer 


Assuming that the source function time ris equivalent to the time t* required for 
stabilization of the ground cloud, the mean standard deviation of the wind azimuth 
angle in the surface mixing layer is 



(2-37) 


For physically reasonable combinations of z , p and t , which are interrelated 

1 Ju 

because of their joint dependence on atmospheric stability in the surface mixing 

layer, the value of cr' {L = 1} falls within the range 
ALi 


^AR 


= 600} 


4 


< or 


AL 


(L = 1} < 


{r, = 600} 


1.5 


(2-38) 


when the reference standard deviation o-' „ is measured near the surface (2 m < z_ 

AR “ R 

< 20 m) over a reference time of 10 minutes. The passage of the rocket vehicle 

through the surface layer introduces trubulence that may persist for 30 minutes or . 

longer, depending on the stability. This vehicle-generated turbulence enhances 

ambient turbulence levels and thus acts to increase the value of 

these reasons, the simple expression 
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(2-S9) 


"11 

2 

is used in the Preprocessor Program to define the mean layer standard deviation of 
the wind azimuth angle. 


The standard deviation of the wind elevation angle is not related to the source 
function time, but is approximately related to the value of {t by the expression 
(Cramer, et al. , 1964) 


o' 


ER 


= 600 seconds} 


3 


(2-40) 


and to the wind profile in the surface mixing layer by 


0. 3 -p 


ER 


R 


unstable 

conditions 


'’e « 



r ( 2 - 41 ) 


; neutral and 
stable conditions 


Therefore, under unstable atmospheric conditions, o' increases with height. Also, 

£ 

while o' decreases with height under neutral and stable conditions, vehicle-induced 
E 

turbulence will persist longer under these conditions and will tend to increase ambi- 
ent turbulent levels. For these reasons, the Preprocessor Program assumes tur- 
bulence is isentropic in the surface mixing layer such that 

o' (t = 600} 

a^j^{L=l} = (2^2) 

2 
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600 


Thus , the only turbulence input required for program operation is cr' {t = 

O 

seconds}, the reference standard deviation of the wind azimuth angle measured over 
a 10-minute period at a height near the surface. At Kennedy Space Center, these 
measurements are available from bi-directional vanes mounted on towers throughout 
the launch complex areas. 

In the NASA/MSFC Multilayer Diffusion Models Program, values must be 
assigned to the lateral a and vertical (3 power-law diffusion coefficients. Experience 
in application of the models to predict dispersion from nearly -instantaneous sources 
has shown that a and j3 normally vary insignificantly from unify. The diffusion 
experiments discussed by Pasquill (1974) indicate that a and for an instantaneous 
source are essentially independent of atmospheric stability. For short travel 
distances, the experimental evidence suggests that o: is greater than or equal to 
unity. Even for travel distances as long as 500 kilometers, a value of unity (Pasquill 
1974, p. 224) appears to provide a reasonably good approximation to the data. Before 
vertical cloud expansion is restricted at the top of the surface mixing layer, the 
diffusion experiments analyzed by Pasquill show that ^ is also equal to unity. At 
larger downwind distances, the experiments indicate a less rapid inc\-ease of cloud 
vertical expansion, which probably reflects the restriction on vertical growth at the 
top of the surface mixing layer. Since the NASA/MSFC Multilayer Diffusion Models 
provide for restricting vertical growth at the top of the surface mixing layer, a value 
for (3 of unity would provide a good estimate of vertical expansion rates. For these 
reasons, the Preprocessor Program sets all power-law diffusion coefficients in the 
surface mixing layer to unity. 

The depth of the surface mixing layer, z (L - 1} for Model 4 and z (K = 1} 

XL i IV 

for Model 3, is not automatically derived from meteorological data available to the 
Preprocessor Program at the present time, but must be specified as an input to the 
program. In the input list to the Preprocessor Program, the depth of the surface 
mixing layer is denoted by An experienced meteorologist can estimate the 
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depth of the surface mixing layer through inspection of the radiosonde data and a 
general knowledge of the mesoscale and ^noptic scale weather patterns at a partic- 
ular site. 
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SECTION 3 


THE NASA/MSFC MULTILAYER DIFFUSION MODEL 

The meteorological structure in the lower troposphere is usually comprised 
of several layers with distinctive wind, temperature and humidity fields. Horizontal 
spatial variations in wind regimes can also occur in the surface layer, usually as a 
consequence of changes in terrain or at land-water interfaces. The models described 
below have been designed to accommodate to these variations of meteorological 
structure in the lower troposphere. The vertical stratification problem in the trop- 
osphere is handled by applying the models to individual layers in which the meteoro- 
logical structure is reasonably homogeneous. Layer boundaries are placed at the. 
points of major discontinuities in the vertical profiles of wind, temperature, and 
humidity. For simplicity, it is assumed that, in general, there is no flux of material 
across layer boundaries due to turbulent mixing. Provision is made, however, for 
the flux of material across layer boundaries as a result of gravitational settling or 
precipitation scavenging and, in Model 4, as a result of breakdown in meteorological 
layer structure. Changes in the meteorological structure of layers, at some 
arbitrary time or downwind distance from the point of release, are accommodated 
by stopping the transport and diffusion processes in the layers affected by the change 
in structure, calculating new sets of initial source and meteorological input parameters 
and re-starting the transport and diffusion process with the new inputs. The provisions 
in the use of Model 4 for permitting changes in initial conditions are also ideally suited 
for calculating concentration and dosage fields in the surface mixing layer when 
portions of the cloud of exhaust products are located at varying spatial positions at 
cloud stabilization. The use of this feature of Model 4 is further explained in Section 
3,4 below. 
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The major changes in the model construct from those appearing in the 
report by Dumbaiild, Bjorklund and Bowers (1973) are that provision has been made 
in the Vertical Terms of Models 3, 4 and 6 for partial reflection of material at the 
ground or water surface. For convenience, the mathematical specifications for all 
the various layer models contained in the NASA/MSFC Multilayer Diffusion Models 
Program are given below. 

3. 1 MODEL 1 

In this layer model, the source extends vertically through the entire layer 
and turbulent mixing is occurring. It is assumed that the vertical distribution of 
toxic material is uniform with height and that the distributions of toxic material 
along the x- and y-layer coordinates are Gaussian. 

3, 1. 1 Dosage Equation for Model 1 


th 

The dosage equation for Model 1 in the K layer is 



In the above expression 

= the source strength in units of mass 

= height of the top of the K layer 

th 

= height of the base of the K layer 

The quantity u^^ in Equation (3-1) is the mean cloud transport 
speed in meters per second in the K layer. In the surface layer (K = 1), the wind 
speed-height profile is defined according to the power-law expression 


Q 


K 


TK 


BK 


32 



(3-2) 


where 


u, 


R 


2j^{K = 1} = 



mean wind Speed measured at the reference height z 

R 

power-law exponent for the wind speed profile in the surface 
layer 



mean wind speed at the top of the surface layer z {K = l} 

TK 

height in the surface layer 


Thus, in the surface layer, the mean cloud transport speed is defined by the ex- 
pression 


Uj^{K = l} = 


C 


R 


(^TK " ^b) 


TK 


dz 


“R 




In layers above the surface layer (K >1), the wind speed-height profile is assumed 
linear and defined by the expression 


u{z^,K>l} 




(3-3) 
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where 


TK 


iTiean wind speed at the top of the layer z 


TK 


BK 


mean wind speed at the base of the layer z 


BK 


, th 

In the K layer (K >1), the mean cloud transport speed is given by the expression 


{k > 1 } 2 


The standard deviation of the crosswind dosage distribution 
is defined by the expression 


yK 
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4.3 


where 


°AK = mean layer standard deviation of the wind azimuth angle in 

radians for the cloud stabilization time r 


In the surface layer (K = 1), 


(“*0(^TK ^r)W” 


(3-5) 
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where 



*^AR 



standard deviation of the wind azimuth angle in radians 
at height z^ and for the cloud stabilization time 


“ar 



standard deviation of the wind azimuth angle in degrees 
at height z^^ and for the reference time period 


m = 


“ATK 


{t, K=l} = 


*^ATK ^o’ 



power-law exponent for the vertical profile of the 
standard deviation of the wind azimuth angle in the 
surface layer 



standard deviation of the wind azimuth angle in degrees 
at the top of the surface layer z for the reference 
time period 


For layers above the surface (K > 1), 
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where 


°ATK 


("aBK t ^ ‘'ABK ( t})/2 


(3-6) 
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standard deviation of the wind azimuth angle in degrees 
at the top of the layer for the reference time period 
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standard deviation of the wind azimuth angle in degrees 
at the base of the layer for the reference time period 


"yK 


downwind distance from the source 
crosswind distance from the axis of the cloud 
crosswind virtual distance 
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when o- {k} > a' { t) x 
yo AK ry 


o 

b 


standard deviation of the lateral source dimension 
in the layer at downwind distance 

X 


distance from the source at which cr {icjis mea- 

Ry 


yo 

sured 

X 

ry 


distance over which rectilinear crosswind expansion 
occurs downwind from an ideal point source 

a 

K 


lateral diffusion coefficient in the layer 


=5 

vertical wind direction shear in the layer 



(®TK " ®Bk)(i8o) 

®TK 


mean wind direction in degrees at the top of the layer 

e 

BK 

= 

mean wind direction in degrees at the base of the layer 

3.1.2 

Concentration Equation for Model 1 


til 

The maximum concentration for Model 1 in the K layer is given by the 
expression 


xK 


(3-7) 


where 


*xK 


standard deviation of the alongwdnd concentration 
distribution in the layer 
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alongwind cloud length for a point source in the layer 
at the distance from the source 
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= ■< 


0.28(|AUjj)(x^) 


{k}< 0 


(3-9) 


’ ^ ®* aI -® 


AUj^{K=1} = 


vertical wind-speed shear in the layer 


AUj^{K>1} = 


u. 


TK 
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BK 


Az 


{K} 


lapse rate of virtual potential temperature in the 
layer 



standard deviation of the alongwind source dimension in 
the layer at the point of cloud stabilization 
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The above equation for L {x } is based on the theoretical and empirical results 
reported by Tyldesley and Wallington (1965) who analyzed ground-level concentration 
measurements made at distances of 5 to 120 kilometers downwind from instantaneous 
line-source releases. 


The maximum centerline concentration for Model 1 in the K layer is given 
by the expression 


^CK ^k} 

The average alongwind concentration is defined as 



( 3 - 10 ) 


^ K pK 


(3-11) 


where 


pK 


th 

cloud passage time in seconds in the K layer 


th 

The time mean alongwind concentration in the K layer is defined by the expression 


^K’ ■* T^ 



(3-12) 


where 

T = time in seconds over which concentration is to be 

A 

averaged 

The time mean alongwind concentration is equivalent to the average alongwind 
concentration when t^^^ equals T^. 
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3.2 


MODEL 2 


Layer Model 2 refers to the same source configuration as Model 1 in 
which the source extends vertically through the entire depth of the layer and the 
distribution of toxic material is uniform with height. In Model 2» however, it is 
assumed that no turbulent mixing is occurring. Consequently, there is no dilution 
of the cloud due to turbulent expansion. The dosage and concentration equations for 
Model 2 are given by Equations (3-1) and (3-7), respectively, with the following 
substitutions: 


"yK ' 


3. 3 MODEL 3 


This layer model differs from Models 1 and 2 in that the vertical extent of 
the source is less than the depth of the layer. The model equation thus contains 
•Vertical expansion terms. 


3.3.1 


the e}q>ression 


Dosage Equation for Model 3 

th 

The dosage equation for Model 3 in the K layer is given by 



(3-15) 
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where 


Q 


K 


source strength or total mass of material in the layer 


R = effective source height or height of the centroid of tiie 

stabilized cloud 

ff „ = standard deviation of the vertical dosage distribution 

zK 

in the layer 


fraction of material reflected at the surface z (1 
for complete reflection and 0 for no reflection) and 
0° is defined to be equal to unity for convenience in 
writing Equation (3-15) 


The remaining terms are the same as those in Equation (3-1) for Model 1. 


The standard deviation of the vertical dosage distribution is 
defined by the expression 
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(3-16) 


where 


EK 


mean standard deviation of the wind elevation angle 
in radians for the layer 


X 


zK 


vertical virtual distance in the layer 

r 


P 


K 


vertical diffusion coefficient in the layer 


X = distance over which rectilinear vertical expansion 

rz 

occurs downwind from an ideal point source 


In the surface layer (K = 1), 
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(3-17) 


where 

°"ER 


standard deviation of the wind elevation angle in degrees 
at the height z 

rC 


q = power-law exponent for the vertical profile of the 

standard deviation of the wind elevation angle in the 
surface layer 
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.^etk = 


standard deviation of the wind elevation angle in degrees 
at the top of the surface layer 


Above the surface layer (K >1), 

°^EK ~ (^ETK ^ '^EBk)(36o) 


where 

a 

ETK 


standard deviation of the wind elevation angle in degrees 
at the top of the layer 


a 

EBK 


standard deviation of the wind elevation angle in degrees 
at the base of the layer 


The vertical virtual distance x „ is given by the expression 
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rz 


standard deviation of the vertical dosage distribution 

at x_ 

Rz 

distance from the source at which o’ (k) is mea- 
th 

sured in the K layer 
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3. 3.2 


Concentration Equation for Model 3 


The concentration equation for Model 3 is the same as that 

for Model 1 which is given by Equation (3-7) in Section 3. 1.2 with D from Equation 

K 

(3-15). Equation (3-10) in Section 3. 1. 2 also gives the maximum centerline concen- 
traion for Model 3. Similarly, average and time mean alongwind concentrations for 
Model 3 are given by Equations (3-11) and 3-12) with from Equation (3-15). 

3. 4 MODEL 4 

Model 4, the layer-breakdown model, may be used to calculate concentra- 
tion and dosage fields resulting from changes in the meteorological layer structure. 
Model 4 may also be used to determine concentration and dosage fields in the surface 
mixing layer downwind from a source in which the source strength varies with height 
in the layer and/or where there are spatial differences in source locations within a 
layer. The application of Model 4 requires the following assumption: 

• The boundary between adjacent layers or sublayers is eliminated 
and the layers are replaced by a single layer L 

• Turbulent mixing is occurring in layer L 

• The material in each of the layers or sublayers is initially uniformly 
distributed in the vertical 

• Reflection occurs at the upper and lower boundaries of layer L 

The selection of Model 4 for layer breal<down calculations or to accommo- 
date vertical source strength variations and/or spatial differences in source loca- 
tions in the surface mixing layer is controlled in the computer program by selection 
of certain options (see Appendix B) available in the input configuration. If no special 
provision is made and Model 4 is specified for use, the program assumes that the 
function of the model is to accommodate to vertical source strength variations or 
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spatial differences in source locations. For example, the surface mixing layer can 

be divided into several sublayers where the source strength, although assumed to be 

th 

vertically uniform in the K sublayer, increases with height in subsequent layers. 

Also, because of wind speed and direction shear in the sublayers, portions of the 
exhaust cloud may be located at different positions in the horizontal plane at the 
time of cloud stabilization. In this case. Model 4 calculates the contribution from each 
sublayer to the composite concentration and dosage fields in the surface mixing layer 
by permitting turbulent mixing across the initial sublayer boundaries. 

If Model 4 is to be used to predict the concentration and dosage fields 

downwind from a change in meteorological structure, the program option 1SKIP(7) 

must be properly set, the input parameter NBK must be initialized, and the meteoro- 

ttl * 

logical parameters for the new L layer and the time t at which layer breakdown 
occurs must be specified (see Appendix B). 

3. 4. 1 Dosage Equation for Model 4 

The dosage equation for Model 4 for the contribution from the 
th 

portion of the cloud in the K layer to the receptor position in the layer L is given 
by the expression 
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erf 
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+ erf 






zLK 


VI <x 


zLK 


where, again for convenience, 0 ° is defined equal to unity. 


The total contribution to a receptor position in layer L is calculated by summing 
the contributions from all K layers. A derivation of the vertical term (terms 
following the summation signs) in Equation (3-18) is presented in Appendix E. In 
the above expression 



u, 


TK 


source strength in units of mass for the source in 
the layer K 

th 

mean cloud transport speed in the L layer 
n 


{(“tk - "bk) (V * “bk)/^| 
n 

XI (“tk - "ek) 

K=1 


(3-19) 


.th 


wind speed at the top of the K layer 



wind speed at the base of the K 


th 


layer 


and n is the number of sublayers in layer L. 
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The crosswind distance from the axis of the cloud to a 

receptor y (defined positive to the right looking downwind) is given by the expression 
Li 

' 'S') '^SK> ®'l • ““ 


where 


xj.y, 


position of the receptor with respect to the origin of 
the reference coordinate system with the y axis 
positive northward and the x axis positive eastward 


^SK’^SK 


coordinates of the cloud centroid in the K layer at 
time t* with respect to the origin of the reference 
coordinate system 


SK 




'SK 


y,-Ogl*cos 6’jj 




coordinates of the source in the layer with 
respect to the origin of the reference coordinate 
system 
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e 


TK 


9 


BK 


wind direction at the top of the K layer in degrees 
wind direction at the base of the layer in degrees 


a , „ in the L 
yLK 


th 


The standard deviation of the crosswind dosage distribution 
layer is defined by the expression 
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(3-22) 


where 

mean layer standard deviation of the wind azimuth 
angle in radians for the effective cloud stabilization 
time T 



In the surface Iyer (L = 1), 






(3-23) 


where 

/T» / = standard deviation of the wind azimuth angle 

ARL^ ^ 

in radians at height z^^ and for time t 
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‘j\RL 


{-"o> = 


m 


L 


m 


L 


^ATL 


{ r, L=l} = 



standard deviation of the wind azimuth angle in 

degrees at height z and for the reference time 

Jn 

period t 
o 

power -law exponent for the vertical profile of the 
standard deviation of the wind azimuth angle in the 
surface layer L = 1 



or.„T ^^o» = standard deviation of the wind azimuth angle in degrees 

ATLi 

at the top of the surface layer z for the reference 

1 Lt 

time period 

For layers above the surface layer (L >1), 






(3-24) 


where 


°ATL ^ ^ (i„ ) (l8o) 
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^ATL ^ ^ ~ standard deviation of the wind azimuth angle in degrees 

at the top of the layer for the reference time period 
T 

O 

" "abl^'^o^ (^) (l^) 

o^.Bj^{Tq} = standard deviation of flie wind azimuth angle in degrees 

at the base of the layer for the reference time 

The wind-direction shear in radians in the layer is 

given by the expression 



where 


TL 


BL 


mean wind direction in degrees at the top of the 
layer z^^ 

mean wind direction in degrees at the base of the layer 
^BL 


th 

The crosswind virtual distance in the L layer due to source 
,th 


(cloud) originating in the K layer is given by the expression 
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where 



X 


.1/“t 


VKL 


ry 


^AL 


{t} 


ry 


ry 






crosswind source dimension in Layer L due to source 

th 

(cloud) originating in the K layer 

1/2 

alongwind standard deviation of the dosage distribution 
th 

in the K layer at time t* 

crosswind standard deviation of the dosage distribution 
th 

in the K layer at time t* 

lateral diffusion coefficient in the layer 

distance over which rectilinear crosswind expansion 
occurs downwind from an ideal point source 



The downwind distance from the point where the change in layer 

th 

structure occurs for the source (cloud) in the K layer to the point where the dosage 

is to be calculated x is given by the expression 
L 

Xl = -(Xj-Xgj^) sme'^-(y. (3-25) 

The standard deviation of the vertical dosage distribution 
a in the L**' layer is defined by the expression 

ZXjK. 
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where 




rz 


zLK 


X 

SL rz 


it) 


(3r26) 


mean standard deviation of the wind elevation angle 
in radians for the layer 

vertical diffusion coefficient in the layer 

distance over which rectilinear vertical expansion 
occurs downwind of an ideal point source 


In the surface layer (L = 1), 


cr ffz {L = 1}’) - (z 

, = ERL Lv TL V V R/ J / 7T \ 

(^TL ^ 


(3-27) 


where 


'ERL 


standard deviation of the wind elevation angle in degrees 
at the reference height z 

X\ 

power-law e^qponent for the vertical profile of the 
standard deviation of the wind elevation angle in the 
surface layer 


log 






standard deviation of the wind elevation angle in 
degrees at the top of the layer z 

TLi 
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Above the surface layer (L > 1) 




ETL EBL 


)fe) 


(3-28) 


where 




standard deviation of the wind elevation angle in 
degrees at the top of the layer z 

TLi 

standard deviation of the wind elevation angle in 
degrees at the base of the layer z 

dLm 


3. 4. 2 Concentration Equation for Model 4 


by the e:q)ression 


The maximum concentration, equation for Model 4 is given 


D u 
LK L 


Pl- ^l) - ^ 


xLK 


(3-29) 


standard deviation of the cloud alongwind concentration 
distribution in the layer 



L = alongwind cloud length of a point source at distance 


where 

^’xLK 
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= 


0 . 28 (auJXl 


“■28(1 *“lI) (’‘l) 


; AQ_ > 0 
Li 


; Au <0; “{l}< 0 
L Az 


; Au < 0 ~ {l} >0 
L Az 


h (3-30) 


AU, 


vertical wind speed shear in the layer 


“tl - “bl 


f w 


lapse rate of potential temperature in the layer 


xKL 


alongwind source dimension in layer L due to source 

til 

(cloud) originating in the layer 


nr 1/2 


The maximum centerline concentration for Model 4 in the L 
layer is given by the expression 


th 


^CLK (^LK* ^LK~®’ \ k ) TERM) 


-1 


= 


L \^°yLK /J 


(3-31) 


The average alongwind concentration at the cloud centerline 


is defined as 


LK 


Dt IT /t T 

LK/ pL 


(3-32) 
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where 


pL 


til 

cloud passage time in seconds in the L layer 

4. 3 cr _ / u- 

xLK L 


The time mean alongwind concentration in the layer for 


averaging time is defined by the expression 


^LK’ ^LK’ T 


u T 
LK I. A 


U >/2 0 - 


(3-33) 


xLK/ 


3. 5 MODEL 5 

This model is used to calculate the amount of material deposited on the 
surface by precipitation scavenging in the K layer. The ground-level deposition 
WD due to precipitation scavenging, for the case in which the vertical distribution 
of toxic material in the layer is uniform with height, is given by the expression 


WD:,{ 


AQ 






(3-34) 




where 


Q 


K 


■'1 

A 


source strength in units of mass for the source in 
layer K 

time precipitation begins 

percent of material removed per unit time 
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The principal assumptions made in deriving the above expression are: 

/ 

• The rate of precipitation is steady over an area that is large com- 
pared to the horizontal dimension of the cloud of toxic material 

• The precipitation originates at a level above the top of the toxic 
cloud so that hydrometeors pass vertically through the entire cloud 

• The time duration of the precipitation is sufficiently long so that the 
entire alongvvind length of the toxic cloud passes over the point x 

Engelmann (see Slade, 1968, pp. 208-221) discusses the general problems of cal- 
culating the amount of material removed by precipitation scavenging and recommends 
values of the coefficient A that may be combined with precipitation rates to obtain 
estimates of total surface deposition. Other useful information may be obtained 
from the proceedings of the 1970 Symposium on Precipitation Scavenging (Engelmann 
andSlinn, 1970), from Pellett (1974) and from Khutson, et al. , (1974). 

I 


When changes in layer structure occur at time t* the contribution to 

til 

ground deposition WD due to precipitation scavenging in the K layer is given by 

XifC 

the expression 


f 1 

WD < X , y , x=0 > = -p= 

liK t 1* L / 


exp 


yLK L 


2-1 s 
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Maximum ground -level deposition at a point (x ,y , z=0) assuming no 

XJ J-i 

previous cloud depletion due to scavenging, can be obtained by setting the second 
exponential term in Equation (3-35) to unity. Total ground deposition is obtained by 
summing the contributions from all layers through which precipitation is falling at 
points on the reference grid coordinate system. The height of the top of the upper- 
most layer throu^i which precipitation is falling Zj.^ must be supplied as an input 
to the computer program. 

The dosage or concentration at a point in space, assuming precipitation 
scavenging occurs, is obtained by multiplying the appropriate dosage or concentra- 
tion equation by the exponential term in Equation (3-34) or (3-35) containing the 
coefficient A . 

3. 6 MODEL 6 

This model is used to calculate the ground deposition due to gravitational 

settling. The basic source configuration is an area source of finite lateral extent 

and unit vertical extent. Other source configurations are treated by summing the 

deposition at the ground resulting from a number of basic sources arranged to 

simulate the desired configuration. The model is essentially a tilted plume model 

in which the effects of wind shear are taken into account. The axis of a particle or 

droplet cloud of a given settling velocity intersects the ground plane at a distance 

from the source and at an angle from the mean surface wind direction that are 

proportional to the total angular wind shear ajid the residence time of the settling 

material in the layers between the source and the ground surface. In any layer, 

the inclination of the cloud axis from the horizontal is given by tan ^ V / u , where 

s 

V is the particle or droplet settling velocity and u is the mean transport wind 

® th 

speed in the layer. In all cases, material released in the K la}'er and dispersed 

upwards by turbulence is assumed to be reflected downwards at the interface of the 
th th 

K and (K + 1) layers. The basic model is used to calculate the ground-level 
deposition pattern for a single value of the settling velocity. The total deposition , 
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pattern is obtained by summing the results for all settling velocities representative 
of the particle or droplet-size distribution of the released material on a reference 
coordinate grid system. 

In the computer program, provision is made for calculating deposition 
from a source which fills the layer in the vertical and for a source in which the 
vertical extent is less than the depth of the layer. These models are described 

below. 

1 

' 3. 6. 1 Gravitational Deposition Model for a Source that Extends 

Vertically Through the Entire Layer 


Ground-level deposition by gravitational settling for a source 

that extends vertically through the entire layer and in which the material is uniformly 

distributed in the vertical is calculated by summing contributions from a number of 

th 

elementary sources in the K layer. Deposition at the surface for a single ele- 
mentary source at height H „ in the layer is given by the expression 
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Ztt a 
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ynK / J 


(3-36) 


where 


f. 


1 


Q 

T 


K 

K 

K 


fraction of particles or droplets with setting velocity V 

s 

source emission rate in layer K (g sec 
source emission time in layer K 

number of elementary sources in layer K for simulating 
a uniform vertical distribution 
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y = lateral distance from the deposition axis of particles 

s 

or droplets with settling velocily V , 

s ‘ 

= R sin 0 
s ^ s 

Rg = radial distance in the horizontal plane from the source to 

a receptor 

angle between the axis of the ground -level deposition 
pattern and the radial connecting source and receptor 
for settling velocity V 

s 

The terms M and N „ are vertical terms that include provision for reflection from 
nK th 

the boundary between the and (K + 1)^^ layers. These terms are defined by the 
expressions 
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s s 
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nK 


mean \vind transport speed in the .layer between H 
and the ground 
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The quantity u is the mean layer wind speed between the 
height and the base of the layer. The following expressions define the mean 
layer wind speeds in the surface layer (K = 1) and the layers above the surface layer 
(K>1): 


V,{K=1> = 



(3-38) 

(Up)(lI^^{K=l) - 







(3-39) 


The mean standard deviation of the wind elevation angle in 

radians in the layer between H and the base of the layer is given by the 

nlv 

expressions 

, /K=n = / 7T \ 
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The vertical diffusion coefficient in the layer between H 


,th 


nK 


and the base of the K layer is given by the terms 




(3-12) 
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The mean lateral diffusion coefficient in the layer between 

H „ and the surface is given by the terms 
nK 



The number of elementary sources f required to simulate 

K. 

a uniformly distributed source m the vertical is given by the expression 

(3-48) 

^ C^TK"/bk)/^K 


where 


vertical separation of elementary sources in the K 
layer 


R a. 


c EH 



R = a constant value depending on the accuracy desired in 

c 

simulating a vertical line source configuration. A 
value of R^ = 0.45 yields deposition estimates that 
are within 10 percent of the true value 
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The computer program for calculating gravitational deposition automatically dis- 

th 

tributes t sources in the K layer with uniform vertical spacing. The height 
K 

H^K the above equations is the height above the ground of each elementary source. 

The angle between the axis of the ground-level deposition 

pattern and the radial connecting source and receptor for settling velocity V is 

s 

defined by the expressions 


•^s " hi ■ I ^ ^ 

“ hi ~ I ^ ®1 


(3^9) 


where 

0^ = mean wind direction at the reference height 

e = angle between north and a line connecting source and receptor 
R ® 



3. 6. 2 Gravitational Deposition Model for a Volume Source 

in the Layer 

th 

For a volume source at height H in the K layer, the 

bK 

ground-level deposition from gravitational settling is given by the expression 


DEP, 


SK 


V 

Vk 


2 

+ f [if ^SK \ 
(sK )J 


(3-50) 


where the subscript SK indicates that the parameters refer to a single souce in the 
til 

K layer. The subset of equations which define the SK subscripted parameters is 

the same as the subset defining the terms in Equation (3-36), except the following 

substitution is made for the term x appearing in Equation (3-37): 

s 
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(3-51) 


= R 


SK 


cos 


^SK ^zSK 


where 

^zSK 


O', 


ESK 


a (SK) 
zo 


the vertical virtual distance for the volume source 



mean standard deviation of the wind elevation angle 

in the layer between H_ and the ground 

SK 

vertical source dimension of the volume source 


In using Equation (3-50), deposition patterns from all values of V„ representative 

oK 

of the particle or droplet size distribution of the volume source are summed on a 
reference coordinate system to obtain the total deposition pattern. 
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SECTION 4 


DESCRIPTION OF THE NASA/MSFC MULTILAYER 
MODEL COMPUTER PROGRAMS 

The NASA/MSFC computer programs described in this technical document 
consist of a Preprocessor Program and the main NASA/MSFC Multilayer Diffusion 
Models Program - Version 5. Both programs are written in FORTRAN IV and are 
designed for execution on a UNIVAC 1108 computer. In addition, the automated 
plotting routines in Version 5 of the program are designed to provide input to the 
Stromberg-Carlson (SC 4020) machine at MSFC. This section briefly describes the 
general characteristics of the two programs. 

4. 1 PREPROCESSOR PROGRAM 

The Preprocessor Program is designed to automatically calculate source 
and meteorological model inputs for use in Models 3 and 4 of the NASA/MSFC 
Multilayer Diffusion Models Program - Version 5. Requisite inputs to the Prepro- 
cessor Program include vertical profiles of meteorological data of the type available 
from radiosonde soundings, surface turbulence data available at KSC from meteor- 
ological towers, and logic information such as the type of launch and vehicle for 
which the calculation is being performed. The program can currently perform 
calculations for the following vehicles: 

• Space Shuttle 

• Titan III 

• Delta-Thor 

• Minuteman II 

Calculations performed by the Preprocessor include: 

• Time-height profile of the rise of hot exhaust products contained in the 
ground cloud 
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• Position in space of the stabilized ground-cloud 

• Source dimensions of the stabilized ground-cloud 

• Distribution of pollutant products within the stabilized ground-cloud 

• Turbulence input parameters 

The cloud-rise models and algorithms used in computing this information are described 
in Section 2. The data decks produced on option by this program include a complete 
card deck for each of the four pollutants HCl, CO, CO 2 and AlgOg for use in Models 
3 and 4 of the NASA/MSFC Multilayer Diffusion Models Program. The program also 
provides the option to produce input data for normal and abnormal launches of all four 
vehicles. A complete description of the program with user' s instructions is presented 
in Appendix A. 

4. 2 NASA/MSFC MULTILAYER DIFFUSION MODELS PROGRAM - VERSION V 

The NASA/MSFC Multilayer Diffusion Models Program is designed to calcu- 
late the following quantities downwind from normal and abnormal launches of rocket 
vehicles: 


• Concentration and dosage patterns 

• Time - mean concentration patterns 

• Average cloud concentration 

• Time of cloud passage 

• Ground -level deposition patterns due to gravitational settling or precip- 
itation scavenging 

Program options include the calculation of concentration, dosage and time-mean 
concentration patterns with partial reflection of material at the surface, with time- 
dependent exponential decay, and/or with depletion due to precipitation scavenging. 
Also, the program is capable of calculating ground-level deposition due to gravitational 
settling with partial reflection at the surface. Provision is also made (in Model 4) 
to account for changes in meteorological structure along the cloud trajectory. A 
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description of the mathematical specifications of the models included in the program 
is given in Section 3 above. 

The multilayer model program is capable of accepting input data from cards 
from the Preprocessor Program or from input card data sets supplied by the user. 

Program output options include; 

• Printing of all data inputs 

• Printing of the results of all model calculations 

• Plotting of maximum centerline concentration, dosage, time-mean 
concentration and deposition versus distance along the cloud trajectory 

• Plotting of concentration, dosage, time-mean concentration and depo- 
sition isopleths 

A simplified block diagram illustrating major features of the program is s 
shown in Figure 4-1. A description of the NASA/MSFC Multilayer Diffusion Models 
Program complete with user' s instructions is contained in Appendix B. 
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Figure 4-1 Simplified block diagram of the computer program for the NASA/ 
MSFC Multilayer Diffusion Models. 
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SECTION 5 


EXAMPLE CALCULATIONS 

Example calculations have been made for both normal and abnormal launches 
of the Space Shuttle vehicle to illustrate the use of the Preprocessor Program in con- 
junction with the NASA/MSFC Multilayer Diffusion Models Program - Version 5. The 
meteorological data used in the example calculations are described in Section 5. 1. 

The results of the calculations for Model 3 are given in Section 5. 2 while those for 
Model 4 are given in Section 5. 3. Computer printout of the example calculations is 
presented in Appendix D. 

5. 1 METEOROLOGICAL INPUT DATA 

Ground-level concentrations and time-mean concentrations of HCl, CO and 
A1„0„ were calculated for normal and abnormal launches of the Space Shuttle using 
meteorological measurements made at Kennedy Space Center on 21 October 1972. 
Surface weather maps show that a cold front approached Florida on 19 October and 
passed the Cape on the morning of 20 October. By 21 October, the cold fi'ont was 
located just south of Florida. Meteorological data from a radiosonde released at 
1115 Z on 21 October are given in Table 5-1. Figure 5-1 shows temperature, wind- 
speed and wind-direction profiles obtained from these data. The temperature profile 
shows temperature decreasing with height to about 1600 meters above the surface 
with the more rapid decrease occurring between 300 and 1400 meters. The wind 
speed increases from 6 meters per second near the surface to 11 meters per 
second at 750 meters, remains constant between 750 and 1432 meters, and then 
decreases with height. As shown in Table 5-1, the relative humidity increases with 
height to 1432 meters and then decreases. From this information, the depth of the 
surface mixing layer H^, which is a required input to the Preprocessor Program, 
was set equal to 1432 meters. The surface density at the time of the radiosonde 
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RADIOSONDE MEASUREMENTS FOR 1115 Z ON 21 OCTOBER 1972 
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launch was 1197. 07 grams per cubic meter. Bi-directional vane measurements 
from towers at KSC indicated that the reference standard deviation of the wind 
azimuth angle 0 -^^^ for a ten-minute period was approximately 9 degrees at a height 
of 18 meters. Other meteorological input parameters required by the Preprocessor 
Program are supplied by the data given in Table 5-1. 

5 . 2 RESULTS OF CALCULATIONS USING MODEL 3 OF THE NASA/MSFC 

MULTILAYER DIFFUSION MODELS PROGRAM 

Table 5-2 contains the results of the Preprocessor Program calculations 
of the Model 3 inputs using Space Shuttle source data and the meteorological inputs 
for 21 October 1972 at KSC. The cloud stabilization height for a normal 
launch was calculated using Equation (2-2) while Equation (2-8) was used for the 
two abnormal launches. Effective heists of the stabilized cloud in the surface 
mixing layer are based on the relationships given by Equation (2-13). The source 
dimensions were calculated from Equations (2-11) and (2-12). Figure 5-2 shows 
the configuration of the exhaust product cloud at the time of stabilization for the 
normal launch. The stippled area in Figure 5-2 represents the dimensions of the 
stabilized cloud in the surface mixing layer used in the Model 3 calculations. The 
position of the stabilized cloud in relation to the launch pad and the time required 
for cloud stabilization are also given in Table 5-2. The source strengths of the 
three exhaust products HCl, CO and AI 2 O 3 were calculated from Equation (2-32). 

The results of the Model 3 calculations of the maximum centerline v and 

c 

ten-minute average x {lO niin} HCl concentrations at the surface for a normal launch, 
single-engine burn and slow-burn on the pad are shown in Figures 5-3 through 5-5. 
Figure 5-3 shows, for a normal launch, a calculated peak maximum centerline HCl 
concentration of 0 . 8 parts per million (ppm) at a distance of about 12. 5 kilometers 
downwind from the launch pad. The profile of centerline (y = 0) ten -minute average 
HCl concentration, which is obtained by calculating the mean concentration over the 
ten-minute period of highest concentrations as the cloud passes a point in space, 
shows a maximum of 0.21 ppm HCl at a distance of 12. 5 kilometers from the 
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TABLE 5-2 


MODEL 3 INPUT PARAMETERS PRODUCED BY THE 
PREPROCESSOR .PROGRAM 


Model 3 

Input Parameters 

Type of Launch 

Normal 

Single-Engine Burn 

Slow Bum 

Cloud Stabilization 

Height, z (m) 
m 

Effective Height, 
H(m) 

a) Cloud Height 


1790 

1038 

1910 

1194 

1772 

1159 

Dimeasion (m) 

cr = O' 

xo yo 

or 

ZO 

b) Source Dimensions 

533 

183 

444 

111 

412 

127 

Time of Cloud 
Stabilization t* 
(sec) 

c) Time of Rise 

447 

400 

452 

Range, R^^ (m) 
Azimuth, A^^ (deg) 

d) Cloud Position from Pad at Stabiliation 

4103 

235 

3646 

230 

4218 

235 

Exhaust Product 

HCl (mg) 

CO (mg) 

AI 2 O 3 (m^ 

e) Source Strength in the Surface Mixing Layer 

8.092 X 109 
1.425 X lO^O 

1.802 X 10 ^® 

8.808 X 10^ 
1. 551 X 10^® 
1.961 X 10^® 

2. 555 X 10^® 
4.499 X 10^® 
5.700 X 10^0 
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FIGURE 5-3. Maximum centerline and ten-minute average x {10 min} 

HCl concentrations at ground-level for the simulated normal 
launch of the Space Shuttle on 21 October 1972 using Model 3. 
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FIGURE 5-4^ . 


3 4 5 6 8 iO 20 

DISTANCE (kilometers) 

Maximum centerline and ten-minute average y{10 min} HCl 
concentrations at ground-level for the simulated single-engine 
burn of the Space Shuttle on 21 October 1972 using Model 3. 





HCl CONCENTRATION (ppm) 


301 



FIGURE 5-5. Maximum centerline and ten-minute average )( (10 min} HCl 

concentrations at ground-level for the simulated slow-burn of 
the Space Shuttle on 21 October 1972 using Model 3, 
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launch pad. The corresponding peak maximum centerline HCl concentration for a 
single-engine launch abort, as shown in Figure 5-4, is 0. 84 ppm at a distance of 
about 13. 5 kilometers from the launch pad. The maximum ten-minute average 
concentration of .21 ppm occurs at 15 kilometers from the pad. Comparison of the 
Model 3 calculations for the slow-burn shown in Figure 5-5 with the corresponding 
results for a normal launch and single-engine burn shows that ground-level concen- 
trations are greater downwind from the launch pad following the slow-burn abort. 
The maximum centerline HCl concentration for the slow-burn case is 2. 7 ppm HCl 
at 14 kilometers from the pad and the maximum ten-minute average concentration 
is 0. 62 ppm HCl at 15 kilometers from the pad. 

Concentration and dosage profiles and isopleths of concentration and dosage 
for HCl calculated using Model 3 can be found in the computer listing of example 
solutions in Appendix D. 2. 

5. 3 RESULTS OF CALCULATIONS USING MODEL 4 OF THE NASA/ MSFC 

MULTILAYER DIFFUSION MODELS PROGRAM 


Table 5-3 presents the results of the Preprocessor Program calculations 
of inputs to Model 4 calculations of HCl, CO and AlgOg concentrations downwind 
from normal and abnormal simulated launches of the Space Shuttle vehicle- on 21 
October 1972 at Kennedy Space Center. • The source strength distribution shown In 
Table 5-3 was calculated following the precedures described in Section 2. 3. 2 
which utilizes Equations (2-27) through (2-34). Source dimensions for Model 4 
were calculated using Equation (2-14) and the source positions at the time of cloud 
stabilization were calculated from Equations (2-15) through (2-26). Figure 5-6 
shows the configuration of the exhaust cloud at the time of stabilization for the 
normal launch of the Space Shuttle used in the Model 4 calculations. The abscissa 
of Figure 5-6 is the range from the lavmch pad without consideration of the offset due 
to differences in the azimuth bearings from the launch pad. The calculated cloud 
st^llization heights and stabilization times for Model 4 are identical to those shown 
in Table 5-2 for Model 3 calculations. 
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PREPROCESSOR PROGRAM CALCULATED SOURCE STRENGTHS, SOURCE 
DIMENSIONS AND SOURCE POSITION FOR MODEL 4 
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FIGURE 5-6. Configuration of the stabilized cloud of exhaust products used with Model 4 in calculations for a 
simulated normal launch of a Space Shuttle vehicle at Kennedy Space Center on 21 October 1972, 




The'results of the Model 4 calculations of the maximum centerline and 

ten-minute average HCl concentrations at the surface downwind from a simulated 

normal launch, single-engine burn and slow bum of the Space Shuttle vehicle are 

respectively shown in Figures 5-7, 5-8 and 5-9. The Model 4 HCl peak centerline 

concentrations x are 0. 6 ppm, . 74 ppm and 2.25 ppm for a normal launch, single- 
c 

engine burn and slow burn on the pad. All of these occur at a distance of approximately 
12. 5 kilometers downwind from the launch pad. Similarly, the maximum ten-minute 
average HCl concentrations are . 17 ppm, 0.18 and 0.53 ppm respectively for a 
normal launch, single-engine burn and slow burn on the pad. All of these occur 
at a distance of approximately 15 kilometers downwind from the launch pad. 

Table 5-4 shows a comparison of the results of the Model 4 calculations 
with those obtained from Model 3. The Model 3 maximum concentrations are 15 
to 30 percent greater than the maximum concentrations calculated using Model 4. 

On the other hand. Model 4 predicts higher concentrations than Model 3 at distances 
close to the launch pad. These differences are a consequence of the two source 
configurations as shown in Figure 5-2 (Model 3) and Figure 5-6 (Model 4). Until 
accurate measurements of the vertical distribution of exhaust products in the stabi- 
lized cloud become available, model calculations of concentrations and dosages 
close to the launch pad are subject to this type of uncertainty. At longer downwind 
distances from the launch pad, as comparison of the Model 3 and Model 4 calculations 
shows, there is close agreement in the predicted concentrations because the effects 
of differing assumptions regarding the source configuration are small. 

Calculated concentration and dosage profiles and isopleths of concentration 
and dosage for HCl obtained by using Model 4 are also presented in the computer 
listing of example solutions in Appendix D. 2. 
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FIGURE 5-8. 
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DISTANCE (kilometers) 


Maximum centerline and ten-minute average .v{10 min} HCl 
concentrations at ground-level for the simulated single-engine 
burn of the Space Shuttle on 21 October 1972 using Model 4. 
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FIGURE 5-9. 
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Maximum centerline and ten-minute average y{10 min} HCl 
concentrations at ground-level for the simulated slow-bum of 
the Space Shuttle on 21 October 1972 using Model 4. 






TABLE 5-4 


COMPARISON OF PEAK MAXIMUM CENTERLINE AND 
TEN-MNUTE AVERAGE HCl CONCENTRATIONS (ppm) PROM 
MODEL 3 AND MODEL 4 CALCULATIONS 



Type of Launch 

Model 

Normal 

Single- 

Engine Burn 

Slow Burn 

Number 




X {10 min) 

Xc 

X {10 min) 

3 

0.80 

0.21 

0.84 

0.21 

2.7 

0.62 

4 

0.60 

0.17 

0.74 

0.18 

2.3 

0.53 
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APPENDIX A 


USER INSTRUCTIONS FOR THE PREPROCESSOR PROGRAM 
FOR USE WITH THE NASA/MSFC MULTILAYER DIFFUSION MODELS 
COMPUTER PROGRAM -VERSION 5 

The Preprocessor Program produces a complete set of data decks for 
input to the NASA/MSFC Multilayer Diffusion Models Computer Program - Version 
5. The Program is specifically designed for use with launches of Space Shuttle, 

Titan III C, Delta-Thor, and Minuteman II vehicles. The data decks produced on 
option by this Program include a complete card deck for each of the four pollutants 
HCl, CO, COo and ALO_ for Models 3 and/or 4 in the NASA/MSFC Multilayer 
Diffusion Models Computer Program - Version 5. 

The Preprocessor Program is written in FORTRAN IV and uses approxi- 
mately 9700^^ locations of computer core storage. The Program requires card 
input, print and punch output and executes in less than 2 seconds on a UNIVAC 1108 
computer. A complete computer listing is given in Appendix C. 1 below. 

A. 1 PROGRAM INPUT PARAMETERS 

The Preprocessor Program requires the input of the following meteorological 
parameters: 

(T - Standard deviation of the wind azimuth angle in degrees 

measured at the first reference height over a 10 minute 
time period. 

p - Ambient air density in grams per cubic meter measured 

at z^, 

z - Height in feet or meters at which the meteorological 

measurements are taken. 


A-1 



e 


Wind direction in degrees at z. 


u - Wind speed in knots or meters per second at z. 

T - Ambient air temperature in degrees Celsius at z. 

P - Ambient air pressure in millibars at z. 

RH - Relative humidity in percent at z. 

The Preprocessor Program also requires control information indicating; 

(1) vehicle type, (2) whether the computer run is for a normal or abnormal launch, 

(3) whether z is in feet or meters, (4) whether u is in knots or meters per second, 

(5) height of the surface mixing layer which must coincide with one of the z inputs 
above, and (6) the model being used and the pollutants for which data decks will be 
produced, 

A. 2 PROGRAM INPUT DATA CARD SEQUENCE 

The first card in the input data deck contains general case titling information 
^d is used for a page heading in the Preprocessor print output and is also punched 
in the output data deck for input to the Multilayer Models Program. The second input 
card contains control information and and p at the surface. An example set of 
data input cards is shown in Figure A-1. The exart^ile problem is described in Sec- 
tion 5 in the main body of this report. 

Data Card 1: 

Columns 1-72 - General data set titling information. If input as 

blanks the program will use the last information 
input 


A-2 


Data Card 2; 


Columns 2-4 
Columns 5-7 


Column 8 
Column 9 

Columns 10-45 

Columns 46-55 
Columns 56-65 
Column 67 

Column 68 


Punch these characters indicating the vehicle type. 
If left blank Titan III C is assumed. 

TTN is the Titan III C vehicle; 

STL is the Space. Shuttle vehicle; 

DTH is the Delta-Thor vehicle; 

MIN is the Minuteman II vehicle. 


Punch YES or leave blank if the run is for a normal 
launch. 

Punch N01 if the run is for an abnormal launch where 
a single engine burns on the launch pad. Not produced 
for the Delta-Thor and Minuteman II vehicles. 

Punch N02 if the run is for an abnormal launch where 
a slow burn on the pad occurs. 

(0 is alphabetic) 

Punch M or leave blank if the heists z are in meters. 
Punch F if the heights are in feet. 

Punch M or leave blank if the wind speed u is in 
meters per second. 

Punch K if the wind speed is in knots. 


Punch the date of the meteorological case or any case 
identification information (optional). 

Punch (see Figure A-1). 

Punch p (see Figure A-1). 

Punch a 1 if output for Model 4 is desired; leave 
blank if not. 

Punch a 1 output for Model 3 is desired; leave blank 
if not. 
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Data Card 2 (Continued); 


Column 69 
Column 70 
Column 71 
Column 72 
Column 73 

t)ata Cards 3 to N-1: 
Columns 1-10 
Columns 11-20 
Columns 21-30 
Columns 31-40 
Columns 41-50 
Columns 51-60 
Column 80 


Punch a 1 if output for HCl is desired; leave blank 
if not. 

Punch a 1 if output for CO is desired; leave blank if 
not. 

Punch a 1 if output for A1„0„ is desired; leave blank 


Punch a 1 if output for CO 2 is desired; leave blank if 
not. Produced only for the Titan III C vehicle. 

Punch a 1 if cloud the trajectory range and azimuth 
bearing are to be calculated; leave blank if not. 


Punch z (see Figure A-1). 
Punch 9 (see Figure A-1). 
Punch u (see Figure A-1). 
Punch T (see Figure A-1). 


Punch P (see Figure A-1). 


Punch RH (see Figure A-1). 


Punch and asterisk (*) if this height z coincides 
with the surface mixing layer height; otherwise, 
leave blank. If not punched on any card, the 
program assumes the last z input to be the mixing 
layer height (See Figure A-1). 



Data Card N: 


The last card in the input data deck must be a blank card. Multiple case 
capability is provided by punching a non-blank character in column 80 of this card. 

If column 80 is non-blank the Program expects another complete case, otherwise the 
program stops. 

The Preprocessor Program assumes the decimal point in all of the abbve 
meteorological parameters to be between the sixth and seventh columns of the 10- 
column field if it is not punched. For example, if the 10-column field 41 to 60 
contained the number A^21033M(Ais a blank), the program would interpret it 
as 1210. 32. To avoid improper alignment of a data value, the decimal should be 
punched. Also, unless it is certain that the cloud rise height will not exceed the 
mixing layer height, it is recommended that at least three heights z be input above 
the height coinciding with the mixing layer height. If the cloud rise height exceeds 
the last height input, the program will stop and ask for data at greater heights. 

A. 3 COMPUTER PROGRAM PUNCH OUTPUT 

The Preprocessor Program will punch a complete card deck for direct input to 
the NASA/MSFC Multilayer Models - Version 5. The first card of each of the ouput 
data decks contains $NAM2 and the last card contains $£ND. The possible data 
decks punched are: 


Inputs for: (1) HCl - Model 4 

(2) CO - Model 4 

(3) COg - Model 4 

(4) AlgOg - Model 4 

(5) HCl - Model 3 

(6) CO - Model 3 

(7) CO2 - Model 3 

(8) AlgOg - Model 3 
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COMPUTER OUTPUT FOR THE EXAMPLE DATA 


A. 4 


The Preprocessor Program first prints constant model parameters for the 
vehicle and prints the input data deck for verification. A complete deck for each 
selected pollutant HCl, CO, CO^ and AlgOg is -printed and punched for Models 3 
and/or 4 of the Multilayer Model. Also, the program lists all calculated parameters 
for Models 3 and/or 4 in a tabularized form for verification of Program calculations. 
A complete output listing for the example case shown in Figure A-1 is given in 
Appendix D. 1. 

A. 5 LINKAGE DIAGRAM OF THE PREPROCESSOR PROGRAM 

Figure A-2 shows the subroutine linkage diagram for the Preprocessor 
Program. Each line terminating at a subroutine name represents a subroutine call. 
The Program also references the FORTRAN library functions ACGS, CCiS, SIN, 
ATAN2 , EXP not shown in Figure A-2. A description of the subroutines shown in 
Figure A-2 is given in Section A. 6. 

A. 6 DESCRIPTION OF PROGRAM SUBROUTINES 

Subroutine VHICLE is the main controlling program in the Preprocessor. 
This routine sets the constants for the particular vehicle desired, inputs all data, 
converts inputs into proper units and controls all calculations and output. 

Subroutine DIM34 calculates the source dimensions for Models 3 and/or 4 
and calculates the effective cloud height (see Section 2, Equations (2-11) through 
(2-14)). 


Subroutine PLUMEl calculates the plume rise for instantaneous sources 
(see Section 2.2. 1). 

Subroutine PLUME2 calculates the plume rise for continuous sources 
(see Section 2. 2. 2). 
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Subroutine DELTXY calculates the cloud trajectory (see Section 2, Equations 
(2-15) through (2-26)). 

Subroutine LEAST calculates the lapse rate of virtual potential temperature 
using a least-squares fit (Equation (2-7)). 

Subroutine CONV converts the output source strength into a form suitable 
for punching variable real numbers in a fixed format. 

Subroutine TURB inserts the turbulence parameters cr^ and in the out- 
put data and calculates t^ = t (Equation (2-4)). 

Subroutine TPZ is a linear interpolation function. 

Subroutine CPHI is used in calculating the virtual potential temperature 
(Equation (2-7)) . 

Subroutine DIST4 calculates the source distribution for Model 4 (see Section 

2. 3. 2). 


Subroutine DIST3 calculates the source distribution for Model 3 (see Section 

2.3.2). 


Subroutine OUTPT prints and punches the output data in a namelist format 
suitable for input to the NASA/MSFC Multilayer Diffusion Models - Version 5. 

Subroutine CONST prints the input parameters and resultant calculations 
in a format suitable for error check and calculation review. 

Subroutine ARGOS is a function reference to ACOS, the UNIVAC 1108 arc- 
cosine function. It is referenced in this way because the IBM 7044 at MSEC uses 
ARCOS for the arc-cosine function and program conversion from the UNIVAC 1108 
to the IBM 7044 is easily completed by removal of this entire subroutine. 
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APPENDIX B 


USER INSTRUCTIONS FOR THE NASA/MSFC MULTILAYER 
DIFFUSION MODELS COMPUTER PROGRAM - VERSION 5 


B. 1 PROGRAM DESCRIPTION 


The NASA/MSFC Multilayer Diffusion Models Program - Version 5 is 
designed to calculate patterns of: 

• Concentration 

• Dosage 

• Time -mean concentration 

• Average cloud concentration 

• Time of cloud passage 

• Ground-level deposition due to precipitation scavenging 

• Ground-level deposition due to gravitational settling 

Program options include the calculation of concentration, dosage and time- 
mean concentration with partial reflection, with time dependent decay, and/or with 
depletion due to precipitation scavenging. Also, the Program is enable of calcu- 
lating ground-level gravitational deposition with partial reflection of material at 
the surface. Other Program options include the printing of all data inputs, the 
printing of all model calculations, the plotting of concentration, dosage and/or time- 
mean concentration on the printer page, and the plotting of concentration, dosage and/ 
or time-mean concentration on the SC4020 plotter at Marshall Space Flight Center. 
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The NASA/MSFC Multilayer Diffusion Models Program - Version 5 is 
written in FORTRAN IV and is designed for use on the UNIVAC 1108 computer at 
Marshall Space Flight Center, Huntsville, Alabama. The Program requires 33000^^^ 
locations of core storage on the UNIVAC 1108 computer. The FORTRAN source 
listing is shown in Appendix C. 2 and a complete list of Program input parameters 
and Program options is given in Section B. 2 and B. 3 below. Also, a linkage diagram 
of the Program subroutines is shown in Section B. 6 and a description of the Program 
subroutines is given in Section B. 7. 

B. 2 PROGRAM INPUT PARAMETERS 

■ This section gives a complete description of all Program input parameters. 
A condensed table of the input parameters is given in Section B. 3. The data input 
format is given in Section B. 4 and an example input coding sheet is given in Section 

B. 5. 


NAMCAS - 72 Hollerith characters of general case identification information. 

This information is printed in addition to the adjusted cloud stabil- 
ization height, range and azimuth bearing and the date and time 
of the run as a title page to the output listing. 

TESTNO - The first 36 Hollerith characters (TESTNO (1) - TESTNO (6)) 
contain the meteorological case information. This information 
is printed in the page heading and plot titles following the words 
"THE METEOROLOGICAL CASE IS". 

Characters 37 through 60 (TESTNO (7) - TESTNO (10)) contain 
the name of the rocket vehicle for use in the page heading, (e, g. 
TITAN IIIQ 

Characters 61 through 72 (TESTNO (11) - TESTNO (12)) contain 
the name of the pollutant only if it is not HCl, CO, C02» or AI 2 O 2 . 
(e. g . NOjj) 
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NPS 


This parameter is used to indicate multiple cases. 


ISKIP 
ISKIP (1) 

ISKIP (2) 


If NPS is set to 0, the Program assumes there is another case 
to follow. 

If NPS is set to 1, the Program assumes this is the last case to 
process. 

Program control option array. 

This option, if set non-zero, indicates patterns of concentration, 
dosage, time-mean concentration, deposition, etc. are to be 
calculated and printed on the polar reference grid system defined 
by XX and YY below. The grid system origin is the vehicle 
launch site and all calculation distances are relative to the origin. 

This option, if set non-zero, is used to calculate maximum center- 
line values of concentration, dosage, time -mean concentration, 
and/or deposition along the cloud trajectory relative to the launch 
site. 

If ISKIP (2) is set equal to 1, the model calculations are printed. 

If ISKIP (2) is set equal to 2, the model calculations are plotted. 

If ISKIP (2) is set equal to 3, the model calculations are both 
printed and plotted. 

The maximum centerline concentration, dosage, time-mean 
concentration and deposition are determined by the use of a spline 
function. At each radial distance (XX) from the origin, the Pro- 
gram determines a curve via the cubic spline that passes through 
each angular (azimuth bearing YY) grid coordinate with the 
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calculated maximum roughly in the midpoint of the curve. The 
Program will then determine the maximum value and output, the 
range and azimuth bearing to that maximum (see Section B. 7, 
Subroutine Spline). 

ISKIP (3) - This option, if set non-zero, is used to calculate isopleths of 

concentration, dosage, time-mean, concentration and/or deposition. 

If ISKIP (3) is set equal to 1, the isopleths are printed. 

If ISI<QP (3) is set equal to 2, the isopleths are plotted. 

If ISKIP (3) is set equal to 3, the isopleths are both printed and 
plotted. 

ISKIP (4) - This option is used only with the calculation of ground-level pre- 

cipitation deposition (Model 5). 

If ISKIP (4) is set non-zero, the maximum possible ground-level 
precipitation deposition is calculated at points downwind from the 
cloud position at the time of the start of precipitation (TIMl), 

These calculations are independent of the elapsed time from 
TIMl to the calculation point. 

If ISKIP (4) is set equal to zero, the calculated precipitation 
deposition at points downwind from the cloud position at time 
TIMl is dependent upon the elapsed time from TIMl to the points. 

ISKIP (5) - This option controls the pollutant name and units printed in the 

page heading and plot legend: 

If ISKIP (5) is set equal to 1, the units of calculated HCl concen- 
tration are in parts per million (ppm) and dosage units are in 
parts per million seconds. 
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If ISKIP (5) is set equal to 2, the uuits of calculated CO concen- 
trations are in parts per million (ppm) and dosage units are in 
parts per million seconds. 

If ISKIP (5) is set equal to 3, the units of calculated CO 2 concen- 
tration are in parts per million (ppm) and dosage units are in 
parts per million seconds. 

If ISIdP (5) is set equal to 4, the units of calculated AI 2 O 2 concen- 

O 

tration are in milligrams per cubic meter (mg/m ) and dosage 
units are in milligram seconds per cubic meter. 

If TESTNO (11) above is non-blank, then ISKIP (5) is used only for 
units selection and the pollutant name is taken from TESTNO (11). 
Also, calculated precipitation deposition (Model 5) and gravitational 
deposition (Model 6) are in units of milligrams per square meter, 

ISKIP (6) - This option is used for printing purposes only and gives the type 

of vehicle launch for which calculations are being made and inserts 
the following in the page heading and plot legend. 

If ISKIP (6) is set equal to 1, a "STATIC FIRE" is assumed. 

If ISKIP (6) is set equal to 0 or 2, a "NORMAL LAUNCH" is 

assumed. 

If ISKIP (6) is set equal to 4, a "SLOW BURN" is assumed. 

If ISICIP (6) is set equal to 5, the program omits this option from 

the page heading and plot legend. 

ISKIP (7) - This option controls the meteorological data used with Model 4. 
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ISKIP (8) 


NSX 


NYS 


If ISKIP (7) is set equal to zero, the Program assumes Model 4 
is being used to determine concentration, dosage, etc. , in a layer 
where the pollutant distribution at cloud st^ilization varies sub- 
stantially with height. The meteorological data used in Model 4 
is automatically determined* from the meteorological inputs 
assigned to the initial layers or sublayers. 

If ISKIP (7) is set equal to 1, the Program assumes Model 4 is 
being used to determine concentration, dosage, etc. , resulting 
from changes in the meteorological layer structure. The meteor- 
ological data used in Model 4 after time TAST (time of layer 
structure change measured from time of cloud stabilization) is 
taken from the input parameters ALPHL through TEMPL (see 
page B-22 below). 

This option, if set non-zero, prints a detailed listing of all 
Program inputs. 

Number of radial distances (range) XX in the polar reference 
grid system. If NXS is set < 0, the default value of 41 is used 
for NXS and the array XX is automatically filled from values 
shown in Table B-3. 

Number of azimuth bearings in the polar reference grid system. 

If NYS is set < 0, (his parameter is automatically calculated 
and the array of azimuth bearing coordinates (YY) is automatically 
filled. The value of NYS includes sufficient points in YY to 
provide a calculation pattern spanning 100 degrees (see Table 
B-3, note 9). 
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NZS 


Total number of initial layer boundaries including the ground 
surface boundary. 


NCI, NDI, NTI - 


NPTS 


NVS 


These parameters each contain two values used in tlie maximum 
centerline calculations under ISKIP (2) and in the calculation of 
isopleths under ISKIP (3). 

The total number of isopleth values is given in the hundreds and 
tens positions of NCI, NDI and/or NTI. If these positions are 
zero, isopleths for the respective quantity (concentration, dosage, 
time-mean concentration and/or deposition) is not calculated. 

The number of critical pollutant levels (air quality standards) to 
be identified in the plots for maximum centerline calculations is 
given in the units position of NCI, NDI and/or NTI. If this 
position is zero, no plot is generated. If set to 9, a plot is 
generated without indicators for critical pollutant levels (air 
quality standards). 

If the units position of NCI, NDI and/or NTI is greater than zero 
and not equal to 9, the critical pollutant levels (standards) must 
be punched as the first values in the arrays Cl, DI and/or TI 
below. 

Number of heights at which calculations are to be made. If NPTS 
is set equal to zero or omitted, NPTS is defaulted to 1 and ZZL (X) 
below is set equal to zero. 

Number of droplet or particle terminal fall velocities used to 
calculate ground-level gravitational deposition from all layers 
except the layer in which a destruct occurs (Model 6 only). 
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NVB 


XX 


YY 


Z 


DELX 


DELY 


Number of droplet or particle terminal fall velocities used to 
calculate ground-level gravitational deposition from the layer in 
which a vehicle destruct occurs (Model 6 only). 

Array of radial distances (range) for the coordinates used in 
calculations on the reference grid system. This array is auto- 
matically filled if NXS = 0, (see NXS above). The last 2 points 
in XX are used only for calculating isopleths; the second to last 
point should equal 1. 2 times the third to the last point and the 
last point should equal 1.5 times the third to the last point. 

Space the XX values uniformly and use as many as the program 
will allow. The user is cautioned to use the default values unless 
another grid is required. 

Array of azimuth bearings for the coordinates used in calculations 
on the reference grid system measured clockyise from zero degrees 
north. This array is automatically filled if NYS = 0 (see NYS above) . 
Space the YY values densely toward the center of the calculation sec- 
tor and use as many as the program will allow. The user is cautioned 
to use the default values unless another grid is required. 

Array of layer boundary heights in ascending order beginning 
with the surface boundary height (the first layer is always the 
surface layer). 

Array of the radial distances (range) from the source location 
(point of cloud stabilization) in each layer to the center of the 
reference grid system (launch site). 

Array of azimuth bearings to the source location (point of cloud 
stabilization) in each layer, measured clockwise from zero 
degrees north. 


B-8 


TABLE B-1 


SOURCE STRENGTH INPUT UNITS 


Model 

Pollutant 

HCL, CO, COg 

m 

1 

1 

2 

2 

1 

2 

3 

1 

2 

4 

1 

2 

5 

2 

2 

6 

2 

2 


Code definition for Table B-1: 


• (3) Q = Q’ 


22.4 


1013.2 


M 273.16 P 
(Concentration output units are parts per million (PPM)) 


• (D Q = Q' 


where 

Q = Source strength in each initial layer 
Q' = Total weight of the material in the layer in milligrams 
T = Surface temperature in degrees Kelvin 
P = Surface pressure in millibars 
M = Molecular weight of the material 
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TABLE B-1 (Continued) 


(Deposition output units for Models 5 and. 6 are milligrams per 
square meter (mg/m^) and concentration output units for Models 
1 through 4 are milligrams per cubic meter (mg/m^)) 
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UBARK 

SIGAK 

SIGEK 

SIGXO 

BIGYO 

SIGZO 

ALPHA 

BETA 

ZRK 


Source strength within each initial layer. The source strength 
input units depend upon the model used and the pollutcant for which 
calculations are being made. Table B-1 gives the appropriate 
input units for each model and pollutant combination. 

Mean wind speed at ZRK followed by the mean wind speed at the 
top of each layer. 

Standard deviation of the wind azimuth angle for reference time 
T at ZRK followed by the standard deviation of the wind azimuth 
angle at the top of each layer. 

Standard deviation of the wind elevation angle at ZRK followed by 
the standard deviation of the wind elevation angle at the top of 
each layer. 

Standard deviation of the alongwind concentration distribution of the 
source in the layer (alongwind source dimension). 

Standard deviation of the crosswind concentration distribution of 
the source in the layer at a downwind distance XLRY from the 
true source (crosswind source dimension). The default value is 
SIGXO. 

Standard deviation of the vertical concentration distribution of 
the source in the layer at a downwind distance XLRZ from the 
true source (vertical source dimension). 

Lateral diffusion coefficient in the layer (default value is 1). 

Vertical diffusion coefficient in the layer (default value is 1). 

Reference height in the surface layer for meteorological measure- 
ments (default value is 2 meters). 
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TEMPK 


TIMAV 

THETAK 

TAUK 

TAUOK 

H 

XRY 

XRZ 

XLRY 


Virtual potential temperature at each layer boundary z. 

This parameter is used only in the calculation of the wind speed 
shear in the layer. If the wind speed shear Is negative and the 
difference between the virtual potential temperature at the top and 
bottom of the layer is also negative, the Program will use the 
absolute value of the speed shear. • If the temperature difference 
is positive or zero, the program will use a wind speed shear of 
zero. If the layer wind speed shear is positive or zero, the 
virtual potential temperature difference is not used. 

Time over which time-mean concentration and average cloud 
concentration are calculated (default value is COO seconds except 
for CO, where it is 300 seconds). 

Mean wind direction at ZRK followed by the mean wind direction 
at the top of each layer. 

Time required for cloud stabilization. 

Reference time for the standard deviations of the wind azimuth 
angle SIGAK (default value is 600 seconds). 

Adjusted cloud stabilization height. 

Distance downwind from a virtual point source over which recti- 
linear expansion in the lateral occurs (default value is 100 
meters). 

Distance downwind from a virtual point source over which recti- 
linear e^ansion in the vertical occurs (default value is 100 meters). 

Reference distance from the true source at which SIGYO is 
measured (default value is zero). 
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XLRZ 

ZZL 

IZMOD 


Reference distance from the true source at which SIGZO is 
measured (default value is zero). 

Vertical calculation heights. This parameter can include any 
heights within the initial layer structure (default value is zero). 

This parameter designates the model number or numbers for’ 
use in each input layer. A brief description of the six Program 
models is given below and a complete mathematical description of 
each model is given in Section 3 of the main body of this report. 
The possible model nmnber combinations input into IZMOD are 
given in Table B-2. 

1 - Model 1, the source extends vertically through the entire 

initial layer and turbulent mixing is occurring. It is 
assumed that the vertical distribution of material is uniform 
with height and the distributions of material along the along- 
\vind and crosswind cloud axes are Gaussian. The digit 1 
is included in the array IZMOD for each layer in which 
Model 1 is to be used. Also, if any digit of IZMOD is 0, 
the Program assumes Model 1 has been designated. 

2 - Model 2 refers to the same source configuration as Model 1 

in that the source extends vertically throu^ the entire 
depth of the layer and the distribution of material is uniform 
with height. In Model 2, however, it is assumed that no 
turbulent mixing is occurring. The digit 2 is included in the 
array IZMOD for each layer in which Model 2 is to be used 
in the calculations (IZMOD = 2, 2, 2, etc.). 

3 - Model 3 differs from Models 1 and 2 in that the vertical 

extent of the source is less than the depth of the layer. The 
model equation thus contains vertical ejq>ansion terms. The 
digit 3 is input to IZMOD for Model 3 (IZMOD = 3). 

~~ B-13 



TABLE B-2 


POSSIBLE INPUT MODEL NUMBER COMBINATIONS 




IZMOD 


1 


PROGRAM ASSUMES CALCULATIONS ARE MADE USING: 


0 


2 


3 


14 


24 


34 


4 


5 


Model X 


Model 1 


Model 2 


Model 3 


Model 1 is used prior to layer transition and Model 4 is 
used after layer transition occurs at time TAST 


Model 2 is used prior to layer transition and Model 4 is 
used after layer transition occurs at time TAST 


Model 3 is used prior to layer transition and Model 4 is 
used after layer transition occurs at time TAST 


Model 4 is used to accomodate to a variation of source 
strength in the layer and layer transition is immediate 
(TAST=0) 


Model 5 is used and the layer structure and source dis- 
tribution is assumed to be that of Model 1 when only the 
digit 5 is given in IZMOD. The digit 5 can be com- 
bined with any of the above digit combinations (145, 45, 
35, etc. ). When a ,5 is combined with any of the above 
digit combinations the Program assumes the layer struc- 
ture and source distribution of that combination are used 
with Model 5. 

Model 6 


The digits under IZMOD can appear in any order. For example, 
14 is the same as 41 and 154 is the same as 415. 
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4 


- Model 4, the layer-transition model, may be used to calculate 
concentration and dosage resulting from changes in the 
meteorological layer structure. Model 4 may also be used 
to calculate concentration and dosage in a layer where the 
pollutant distribution at cloud stabilization varies substan- 
tially with height. 

The application of Model 4 requires the following assumptions: 

• The boundaries between adjacent initial layers or sublayers is 
eliminated (at time TAST) and the layers are replaced by a single 
layer 

• Turbulent mixing is occurring in the resultant single layer 

• The material in each of the initial layers or sublayers is (before 
time TAST) uniformly distributed in the vertical 

• Reflection occurs at the upper and lower boundaries of the resultant 
single layer 

If the parameters TAST and ISKIP (7) are both set to zero (or omitted from 
the inputs) and Model 4 is specified for use, the program assiunes the function of the 
model is to accommodate variations in the pollutant distribution with height in the 
layer at cloud stabilization. For example, the surface mixing layer can be initially 
divided into several sublayers where the source strength, although assumed to be 
vertically imiform in each sublayer, varies from layer to layer. In this case the 
initial layers are immediately reduced to a single layer and Model 4 calculates the 
contribution from each of the initial sublayers to the composite concentration and 
dosage field by permitting turbulent mixing across the initial layer boundaries. 

IZMOD would contain the digit 4 for each of the respective initial sublayers that 
comprise the resultant single layer. 
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If Model 4 is to be used to predict the concentration and dosage fields 
downwind from a change in meteorological structure, the meteorological parameters 
of the new resultant layer or layers must be specified. Also, the parameter 
ISKIP (7) must be set equal to 1 and the parameter TAST set equal to the time (after 
cloud stabilization) at which the layer transition (meteorological structure change) 
occurs. Each of the initial sublayers that are to be included in a single layer after 
layer transition are specified by including the digit 4 in the array IZMOD. For 
example, assume layers 1 through 4 are to be reduced to a single layer after layer 
transition and layers 5 and 6 are also reduced to a single layer. The first four 
values of IZMOD would include a 4, but they would also include the number of the 
model to be used prior to layer transition (14, 24 or 34). The values of IZMOD (5) 
and (6) for layers 5 and 6 would include a 9 and 4, respectively. The 9 is a special 
flag to separate the resultant 2 layers after layer transition. Also, these last two 
values would include the model number to be used prior to layer transition (14, 

24 or 34). If Model 1 was to be used with 4 in the above example the IZMOD inputs 
would be coded as IZMOD = 4, 4, 4, 4, 9, 4 (or IZMOD = 4*4, 9, 4, or IZMOD = 
4*14, 19, 14, etc.). 

5 - Model 5 is used to calculate the amount of material on 

the surface by precipitation scavenging. The digit 5 must 
be included in the array IZMOD for each initial sublayer 
through which precipitation is occurring. Model 5 uses the 
layer structure and source distribution defined by any one 
of Models 1 through 4. Thus, the array IZMOD must 
include the appropriate model number for each layer that 
describes the layer structure and source distribution. For 
example, assume that Model 4 is being used to accommodate 
to variations in the pollutant distribution with height in the 
surface mixing layer at cloud stabilization and that the 
surface mixing layer has been divided into 6 initial sub- 
layers in which the distribution of material can be 
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considered uniform. Also, assume that precipitation is 
occurring through all 6 layers. The array IZMOD would 
then contain six values equal to 45 for each layer 1 throu^ 

6 (IZMOD = 6*45). 

Model 6 is used to calculate the surface deposition due to 
gravitational settling. The basic source configuration is 
a volume source of finite lateral extent and unit vertical 
extent. Other source configurations are treated by summing 
the deposition at the ground resulting from a number of 
basic sources arranged to simulate the desired configuration. 
The model is essentially a tilted plume model in which the 
effects of wind shear are taken into account. The axis of a 
particle or droplet cloud of a given settling velocity inter- 
sects the ground plane at a distance from the source and at 
an angle from the mean surface wind direction that are 
proportional to the total angular wind shear and the residence 
time of the settling material in the layers between the source 
and the ground surface. In any layer, the inclination of the 

cloud axis from the horizontal is given by tan ^ V / u, 

s 

where V is the particle or droplet settling velocity and u 
s 

is the mean transport wind speed in the layer. In all cases, 

til 

material released in the K layer and dispersed upwards by 

turbulence is assumed to be reflected downward at the inter- 
th til 

face of the K and (K + 1) layers. The basic model is used 
to calculate the ground-level deposition pattern for a single 
value of the settling velocity. The total deposition pattern 
is obtained by summing the results for all settling velocities 
representative of the particle or droplet-size distribution of 
the released material on a reference coordinate grid system. 
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Only IZMOD (1) need be set equal to 6 as no other model 
can be executed in the same case. 


DECAY 

i 

ZLIM 


BLAMDA 


TIMl 

Cl, DI and TI - 


Coefficient of time-dependent decay. If DECAY is set > 0, then 
concentration, dosage, time-mean concentration, etc. , are 
calculated with decay (Does not effect Model 5 or Model 6). 

This parameter is the maximum height through which precipitation 
can occur. If Model 5 is selected, ZLIM is automatically deter- 
mined from IZMOD. If concentration, dosage, etc. , are being 
calculated with precipitation occurring (BLAMDA > 0. 0), 

ZLIM is equal to the upper boundary of the uppermost layer in 
which precipitation occurs (ZLIM is defaulted to Z(NZS)). 

Precipitation scavenging (washout) coefficient. If Model 5 is 
selected, this parameter must be greater than 0. Also, if Model 
1, 2, 3 or 4 is selected with BLAMDA > 0 and without Model 5, 
the Program assumes concentration, dosage, etc. , are to be 
calculated with precipitation occurring. 

Time of start of precipitation measured from the time of cloud 
stabilization. 

Arrays of concentration, dosage and time-mean concentration 
values respectively for which isopleths are calculated. There 
can be two groups of data in each of these arrays, where both 
of the groups are arranged in descending order. The values in 
the first group are critical pollutant levels (air quality standards). 
The number of values in this group is given in the units position 
of the parameters NCI, NDI and NTI respectively. The second 
group of values includes all other isopleth levels desired. The 
total nimiber of values in Cl, DI and TI is given in the hundreds 
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TAST 


GAMMAP 


VS 


PERC 


ACCUR 


and tens positions of NCI, NDI and NTI respectively. If pre- 
cipitation, deposition or gravitational deposition is being calculated, 
the array DI is used for these quantities. 

Time of layer structure change (Model '4) measured from the time 
of cloud stabilization. 

This parameter is 1 minus the fraction of material reflected at 
the surface (partial reflection). If this parameter is set to 0, the 
Program assumes complete reflection; if set equal to .4, 60 per- 
cent (. 6) reflection is assumed; and, if set equal to 1, no reflection 
is assumed. If Model 6 is selected and partial reflection is 
desired, the array GAMMAP must have a value for each particle 
settling velocity category. For all otlier models, only GAMMAP (1) 
need be set. 

Droplet or particle terminal fall velocity distribution used in all 
layers except a layer in which a vehicle destruct occurs (Model 6 
only). 

Frequency of occurrence of each velocity category VS (Model 6 
only). 

Accuracy constant for the line source simulation used in Model 6. 

A value of 0.45 ensures that the calculated ground deposition is 
within 10 percent of the deposition expected from a vertical line 
source. If ACCUR is set to 0. 32, the calculated deposition is 
within 5 percent of that expected from a vertical line source. 
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VB 


PERCB 


SCL 


XMAXEN 


YMAXIN 


XSIZE 


YSIZE 


RASTIN 


Droplet or particle terminal fall velocity distribution used in the 
layer in which a vehicle destruct occurs. The layer must be the 
top layer (Model 6 only). 

Frequency of occurrence of each velocity category VB (Model 6 
only). 

Map scale factor in inches for isopleth plots. If the map scale 
factor is 1 inch = 24000. inches, SCL would be input as 24000. 

If set to zero, the Program will scale the isopleths within the 
boundaries defined by XSIZE and YSIZE below. 

Maximum alongwind distance from the launch site in meters for 
isopleth plots. If set to zero, the Program will use XX(NXS-2) 
as the maximum distance. 

Maximum crosswind distance for isopleth plots in meters. If 
set to zero, the Program will calculate YMAXIN. 

The number of raster counts on the SC4020 in the X or east -west 
horizontal plot axis for isopleths. If set to zero, the Program 
will use 937. 

The number of raster counts on the SC4020 in the Y or north- 
south vertical plot axis for isopleths. If set to zero, the Program 
will use 899. 

The number of raster counts per inch on the SC4020 for isopleth 
plots. If input as zero, the Program uses 163.2. 
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XCIZE 


YCIZE 


XMAXJN 


YMAXJN 


ISW 


JSW 


The number of raster counts on the SC4020 on the X or alongwind 
horizontal axis for maximum centerline plots. If set to zero, the 
Program uses 937. 

The number of raster counts on the SC4020 on the vertical axis 
for maximum centerline plots. If set to zero, the Program uses 
899. 

Maximum alongwind distance in meters from the launch site for 
maximum centerline plots. If set to zero, the Program uses 
XX(NXS-2). 

Maximum number of log cycles for the vertical axis of the 
maximum centerline plots if ISW below equals 0 or 2. Maximum 
value of the vertical axis if ISW below equals 1. If set to zero, 
the Program determines YMAXJN. 

Maximum centerline plotting flag. If ISW is set to 0 or 2, the 
Program plots maximum centerline versus distance on a log-log 
plot. If set to 1, the plot is linear on both axes. 

- Isopleth plot switch. If JSW is set equal to 0, the Program will 
fit a cubic spline function to the discrete isopleth points and plot 
a smooth curve through the points. If JSW is set equal to 1, the 
Program will not use the spline function but will plot straight 
lines between adjacent calculated isopleth points. This option 
has been included because the spline function sometimes fails to 
fit the data points when the isopleths are sharply curved. These 
cases are recognized by a high frequency oscillation along the 
plotted curve and can be corrected by smoothing the curve by 
hand or replotting with JSW set equal to 1. 
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The layer step change (transition) parameters below are used only if, ISKIP (7) 
equals 1 and Model 4 has been selected. These parameters are used only when Model 
4 is being used to predict the concentration and dosage downwind from a change in 
meteorological structure (see IZMOD, Model 4 above). 

ALPHL - Lateral diffusion coefficient in each new layer (Default value is 1) . 

BETL - Vertical diffusion coefficient in each new layer (Default value is 1). 

TAUL - Time required for cloud stabilization in the new layers. 

TAUOL - Reference time for the standard deviation of the wind azimuth 

angle SIGAL in the new layers (Default value is 600). 

ZRL - Reference height in the surface layer for meteorological measure- 

ments. This must be set only if the new bottom layer includes the 
initial surface layer (Default value is 2). 

UBARL - Mean \vind speed at the bottom and top boundaries of each new 

layer. These values are input in ascending order of new layers 
with the value at the top boundary preceded by the bottom. If the 
new bottom layer contains the initial surface layer, UBARL at ZRL 
should be input as the bottom value of this layer. 

SIGAL - Standard deviation of the wind azimuth angle for reference time 

T at the bottom and top boundaries of each new layer. If the 
oL 

new bottom layer contains the initial surface layer, SIGAL at 
ZRL should be input as the bottom value of this layer. 

THETAL - Mean wind direction at the bottom and top boundaries of each new 
layer. If the new bottom layer contains the initial surface layer, 
THETAL at ZRL should be input as the bottom value of this layer. 

TEMPL - Virtual potential temperature at the bottom and top boundaries of 

each new layer. 
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B.3 


CONDENSED TABLE OF INPUT PARAMETERS 


The data input parameters required for the computer Program are given in 
condensed form in Table B-3. The information categories in the table are defined 
as follows: 


NAMELIST 

FORTRAN 

MODEL 

UNITS 

LIMITS 

VALUE 

ARRAY SIZE 


Name of the FORTRAN NAMELIST list to which the variables 
belong. 

Fortran symbolic notation defining the program input. 
Mathematical notation corresponding to the FORTRAN notation. 
Dimensional units of the input parameters. 

Numerical limits on input values. 

Default value should the parameter have a value of 0. 

Maximum number of core locations for the input parameter. 


B.4 DATA INPUT FORMAT 

This Program uses the FORTRAN NAMELIST method to input data. Input 
data must be in a specific form in order to be read using a NAMELIST list. The 
first character in each card to be read must be blank. The first card in the NAME- 
LIST list contains the NAMELIST name NAM2 preceded by the character $ or &. 

The last card in the NAMELIST list contains $END (&END) to terminate the list. 

The form of the remaining data items in the list may be: 

a. VaHdble Name = Constant - The vaHdbU name may be a sub- 
scripted array name or a single variable name. Subscripts must be integer con- 
stants. The constant may be integer, real or Hollerith (nH alphanumeric char- 
acters) data. 
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TABLE OF INPUT PARAMETERS 
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TABLE OF INPUT PARAMETERS 
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TABLE OF INPUT PARAMETERS 
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ZLIM is automatically calculated if IZM(^D contains a 5 (Model 5). 


b. Awtcu/ Name. = Set oi Con&tatvU (separated by commas) - The 0 M.ay 
name, is not subscripted. The iet comtanti consists of constants of the type integer 
or real. The number of constants must be less than or equal to the array size. 

Successive occurrences of the same constant can be represented in the form k* 
aomtant. 


The sequence of the input data parameters within the list is not significant. 

A more detailed e.xplanation of the FORTRAN NAMELIST can be found in most 
FORTRAN language manuals. Section B. 5 shows an example input data coding 
sheet. All Program input parameters are set to zero prior to input of the first 
case. Parameters that are not used or have default values need not appear in the 
input deck. When multiple cases are stacked, all parameters retain their values 
from the last case and are changed only by input. 

B. 5 SAMPLE PROBLEM INPUT DATA 

The sample problem is for the hypothetical launch of a Space Shuttle vehicle. 
The meteorological parameters were measured at Kennedy Space Center for the 
21 October 1972 case described in Section 5. The parameter input sheet is shown 
in Figure B-1 and the Program output is shown in Appendix C.2. The data input 
parameters shown in Figure B-1 were automatically calculated by the Preprocessor 
Program and are shown in the Preprocessor output in Appendix D. 1. Only the HCl 
normal launch cases with Model 3 and Model 4 are given here. The first namelist 
deck shown in Figure B-1 is for Model 4 and the second is for Model 3. 

B. 6 LINKAGE DIAGRAM OF THE NASA/MSFC MULTILAYER DIFFUSION 
MODELS PROGRAM - VERSION 5 

Figure B-2 shows the subroutine linkage diagram for the NASA/MSFC 
Multilayer Diffusion Models Program - Version 5. Each line terminating at a 
subroutine name represents a subroutine call. The asterisk indicates subroutines 
in the SC4020 plot packs^e used only in the UNIVAC 1108 program copy at Marshall 
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(Main Program) 


* 







MSFLD APRNTV* LABLV* 

FIBURE B-2. Linkage diagram of the subroutines in the NASA/MSFC Multilayer Diffusion Models Computer Program - 
Version 5. 


Space Flight Center, Huntsville, Alabama. The program also references the 
FORTRAN library functions SIN, COS, ATAN2, EXP, SQRT, ALOG and ALOGIO 
not shown in Figure B-2. A description of the routines shown in Figure B-2 is 
given in Section B. 8. 

B. 7 DESCRIPTION OF PROGRAM SUBROUTINES 

Subroutine MODEL is the main controlling program in the NASA/MSFC 
Multilayer Diffusion Models Program - Version 5. This routine controls the calcu- 
lations for Models 1 through 4 and passes control to subroutines that calculate 
Models 5 and 6. 

Subroutine BREAK completes all calculations for Models 1 through 4 
(Section 3. 1 to 3,4). 

Subroutine READER is the data input routine for the Program. All input 
default values are set in this routine and most meteorological layer parameters are 
calculated in this routine. 

Subroutine DEPOS controls all calculations for Model 6 (Section 3. 6). 

Subroutine SGP completes all calculations for Model 6 (Section 3. 6). 

Subroutine WASHT completes all calculations for Model 5 (Section 3.5). 

Subroutine TESTR is used to determine special indices used in calculations 
with MODEL 4. This routine also calculates the downwind distance at which the 
layer transition for Model 4 begins. 

Subroutine IDENT is used to initialize the SC4020 plotting routines used at 
Marshall Space Flight Center, 
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Subroutine ENDJOB is used to close out the SC4020 plotting routines at 
the end of a job. 

Subroutine RB8 is used in the calculation of the wind speed and the standard 
deviation of the elevation and the azimuth wind angles in the surface layer (q, m, p, 
ql, ml and pi in Sections 3. 1 to 3.4). 

Subroutine ERTRAN is a time and date routine. This routine is a UNIVAC 
system utility routine that returns the data and time. The following statement is 
used to retrieve the 6 character date (MMDDYY) in the first variable and the 6 
character time (HHMMSS) in the second variable: 

CALL ERTRAN (9, NTFB(l), NTFB(2)) 

Subroutine ISO evaluates the normal error function used in the calculation 
of Model 4 (Equation 3-18, Section 3.4). 

Subroutine COORD calculates the downwind and crosswind distances of a 
grid system calculation point relative to the alongwind cloud axis and the layer source 
location. 


Subroutine RBll is used in the calculation of the wind speed and the standard 
deviation of the elevation and azimuth wind angle in the surface layer. (Equations 
(3-2), (3-5), (3-17), (3-19), (3-23), and (3-27), Section 3. 1 to 3.4). 

Subroutine SIGMA calculates the crosswind, alongwind and vertical standard 
deviations of the dosage distribution for Models 1 through 5. 

Subroutine GENPRT controls the printing and/or plotting of all model 
calculations. 

Subrouting LABELS generates the page heading information for all calcula- 
tions. 
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Subroutine ISSOPT plots the isopleths of concentration, dosage, deposition, 
and/or time-mean concentration on the SC4020 plotter. 

Subroutine LLPLOT plots maximum concentration, dosage, deposition, and/ 
or time-mean concentration as a function of distance in the form of a printer plot on 
the UNIVAC 1108. 

Subroutine PRTTTL prints all page heading information. 

Subroutine NMBRS converts a real number to a sequence of BCD characters 
for insertion into the page heading and for plotting numeric values. 

Subroutine LSSOPT plots maximum centerline concentration, dosage, 
deposition, and/or time-mean concentration as a function of distance on the SC4020 
plotter. 

Subroutine PACKS removes extra blank characters from the page heading 
and packs it into a form suitable for printing. 

Subroutine MSFLD extracts a bit string (byte) from one word and stores it 
jinto another. 

Subroutine HEDING sets up key flags that tell the LABELS subroutine 
whether the label is concentration, dosage, label or other. 

Subroutine INPTS is used to set indices for reading or writing records 
from mass storage. 

Subroutine INTOUT writes or reads records from mass storage. 

Subroutine VRTCLE determines the vertical axis label for maximum 
centerline plots. 
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Subroutine FSTPLT plots information identifying the plots for a set of input 
data decks. 

Subroutine ILAXES draws and labels the plot grids for inaxinium cento rllne 
and isopleth plots. 

Subroutine ILPLOT draws and labels the maximum centerline curve and the 
isopleth curve. 

Subroutine MAXMIN determines the maximum and minimum plot values for 
maximum centerline plots if they are not input. 

Subroutine BOUNDS determines whether the plot curve leaves or enters the 
plot boundaries and interpolates for the intersection of the curve and boundary. 


Subroutine SPLINE calculates the coefficients of a third-degree natural 
spline function used to calculate maximum centerline values and to plot maximum 
centerline and isopleth curves. The function derives its name from the fact that such 
a curve approximates the behavior of a mechanical spline used by draftsman to draw 
a smooth curve through a set of points. The curve has the properties of having 
continuous first and second derivatives, and of being the "smoothest" curve through 
a given set of points in the sense that it minimizes 



cr = 

K 

[sNx)]^dx 

where 


s = 

dosage, concentration, etc. 


X = 

azimuth bearing or distance. 


k 

derivative = 1, 2 
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The function s(x) is the unique curve with these properties and is a piece wise 
function given by a polynomial of maximum degree of 3 in each of the intervals 
[x., X. + l), in general by a different polynomial in each such interval. For a more 
detailed explanation of the spline function, the reader is referred to the book by 
T. N. E, Greville, Theory and Applications of Spline Functions . 

Subroutine LINE2V draws a line between 2 points (SC4020). 

Subroutine PRINTV prints alphanumeric information horizontally (SC4020). 

Subroutine FRANIE V advances the camera frame (SC4020) . 

Subroutine LABLV prints a number (SC4020). 

Subroutine APRNTV prints alphanumeric information vertically (SC4020). 

Subroutine SETMIV sets the plot margins (SC4020). 
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APPENDIX C 


FORTRAN SOURCE LKTINGS FOR THE PREPROCESSOR 
PROGRAM AND THE NASA/MSFC MULTILAYER 
DIFFUSION MODELS PROGRAM - 
VERSION 5 


C. 1 FORTRAN SOURCE LISTING FOR THE PREPROCESSOR PROGRAM 


This section contains the complete FORTRAN source listing of the Pre- 
processor Program. 
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subroutine LEAST(Z*TVrQPD2»K»ISWrZP»TVP) LSTOOiOO 
DIMENSION Z(1)»TVU) UST00200 
IF (K .LE. 1) GO TO 30 LST00300 
L 5 K LST00400 
TVti = 0«0 LST00500 
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FUNCTION aRCOS(A) ACSOOlOO 
ARCOS = ACOS(A) ACS00200 
RETURN ACS00300 
ENu ACS00400 
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C. 2 FORTBAN SOURCE LISTING FOR THE NASA/MSFC MULTILAYER DIFFU- 

SION MODELS PROGRAM - VERSION 5 


This section contains the complete FORTRAN source listing of the NASA/ 
MSFC Multilayer Models Program - Version 5. 
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IF (KSW(l) .L£. 0) GO TO 10 6PT04100 
42* m = 19 GPT04200 
43* GO TO 20 GPT04300 
44* 10 IF (KSW(2) ,UE. 0) GO TO 30 CPT04400 
45* ITI = 31 6PT04500 
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205+ C PLOT maximum CENTERLINE CALCULATIONS GPT20500 
206* IK (ISKIP(2) ♦EQ. 1) 00 TO 480 GPT20600 
207+ HT = 0,0 GPT20700 
208+ XPL = 0,0 GPT20800 
209* YPL = 0,0 GPT20900 
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205 * 48 continue RDR 20500 
ZOb* NTAL = 2 *NBK RDR 20600 
207 * DO So I= 1 »NTAL ROR 20700 
208 * C CHECK MINIMUM VALUES RDR 20800 
209 * IF (SIGAU(I) .LT. , 5 ) SIGAL(I) = .5 RDR 20900 
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164» ALPHNK(I2) = ALPHA (N) SGP16400 
165* 95 continue SGP16500 
166* RETURN SGP16600 
167* END SGP16700 
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SUbRoUTlNE IS0(MR»MT) ^ iSOOOlOO 

common /PaRAMT/ TESTN0{12)» ISKIPdS) »NXS»NYS»N2S»N0I»NCI» IS000200 

ilMtiK»NPTS»NVS»NVB»XX(4l) »YY(4l) #2(16) »DELX(15) »DELY(15) »Q(15) » IS000300 

2UBARK(16) »SIGAK(16) »SIGEK(16) »SIGX0(15) »SIGY0(15) »SIGZ0(15) » IS000400 

3ALPHA(20) »BETA(20) »2RK»TIMAV»THETAK(16) »TAUK»TAU0K»H»XRY»XR2f IS000500 
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4U UO DO 120 I=l»3 , VRT04100 
42* , 120 ISTOk(NCV+I) = UMT(I3+I) VRT04200 
43* NCV = NCV+3 VRT, 04300 
44* 130 CALL PACKS (ISTOR»NCV) VRT04400 
45» RETURN VRT04500 
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40 CONTINUE SPL03400 

RETURN ' SPL03500 

END \ SPL03600 


SUBROUTINE ISSOPT (X, Y» NX, NY» FI ♦ LEGENOr NCHARt DP,NYSi XX,DR, I I , IREC, YISSOOlOO 
IT, I5W3, KLINE, NCVfJM, DECAY, lambda) ISS00200 

COMMON /PLTISO/ SCL , XMAXIN, YMAXIN, XSIZE, YSIZErRASTTN, JSrt ISS00300 

COMMON /BNOS/ XrIT , XLFT , YBOT , YTOP , XPL , YPL ISS00400 
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APPENDIX D 


PROGRAM OUTPUT FOR THE PREPROCESSOR PROGRAM 
AND THE NASA/MSFC MULTILAYER DIFFUSION 
MODELS PROGRAM — VERSION 5 

D. 1 PREPROCESSOR PROGRAM EXAMPLE OUTPUT 

The Preprocessor Program produces both printed output and punched data- 
card output. The Program first prints a detailed list of all inputs including the con- 
stants pertaining to the rocket vehicle for which calculations are being made. The 
first page of the example output listing below shows the constants used for the example 
problem described in Section 5 of the main text for a simulated normal launch of a 
Space Shuttle vehicle at Kennedy Space Center on 21 October 1972. The Program 
next prints the FORTRAN NAMELIST data which are subsequently punched on cards 
for Model 4 calculations using the Multilayer Diffusion Models Program. In this 
example, the first printed list and punched card deck are for calculation of HCl con- 
centrations and dosages; the second list and punched card deck are for CO concentra- 
tions and dosages; the third list and punched card deck are for Al O concentrations 
and dosages. A printed list and punched card deck for the calculation of CO con- 

Cl 

centrations and dosages are not produced for the Space Shuttle vehicle. The Pro- 
gram output for the Model 4 calculations ends with a summary table of the Program 
calculations. The Program then produces similar printed listings and sets of punched 
card decks in the same format for use in the Model 3 calculations. 


D-1 



eg 

n 


biOOOOOOOOOOOOOOOOOOOO 

UOOOOOOOOOOOOOOOOOOOO 

QCOOOOOOOOOOOOOOOOOOOO 

m**»*»»*»************ 

Q^^(vp.^a}uiNo^^*o^r^r*ao'iO(no^^^ 

^lnlnlnvng>g)^»Ncoc^(^c^(^c^r^tnm^^ln 


U|«-00000000000000000000 

aciftoo>oo'o^o*^*-4o\no'K>f^o^*H^O'0^ 

^j;oo%or«-^o^oh*«oocgo*^^oK>iOfO%oio 

I/) c\jOTQC'incoor-oo\o^mcoo\0<\icoio«H 
uiCL.cgoo^a)g}minro«^oo^cog)U)tnio*oocog3 
a: ooo> 0 ' 0 ^ 0 'a^o^o^ 0 'to«oco<cf( 0 <x)<o(t>f^i^ 

CL ^ -« 


:«i^oooooooooooo.ooooooo( 

>0OOOOOOOOOOOOOOOOOOO< 

^^)<g•HO^Oloco<^'0^•ro^o^•tf>roo'rolna-l« 

c ••••••••••••♦•••••• 

:uJcgcgcgcgO‘cor^%o^»oio<goooo*ocg-<( 


txJ 

UJ JO 
uj or* 


X <g bi < 
j <o o »- 
:qco o < ; 

ae . 

:i/i • uji 
4 liJ ^ 1 

w O O 1 
toe ^ O 03 1 
: otneg o < < 

j j -4 • 

- < • u4 

j oo: X I 

: 'all in >“ : 

C X h- »-• _l ( 

: will u u. j 
t xq: X o < 

Ll- < W 
3 X <t X ! 
JQ. < \J o : 

X U.X w I 
<13 O < W 
COW 31 
2 Win an 
J WUJ-t-i w 
i< X ui < r • 
> uix X w w : 
-• X w X mj 
: < o 111 w I 

J JXW W Q I 

► < W Ul < 
W2 wX J 

I OuJV^ <( 
w lU W 2 ' 
1 I X 2 O . 
I 4/1 w : 

> < w w < 
: w X 1 3 0 1 
C <X -J < . 

U1 <CCL O X ( 

c X o oo. u. : 


o ui< w 

X J'X UJ 

< w X UJ 
O w o Ll 

W 2W 

< 


w ^ eg 
^ to 
o ^ 

1 2 xin w 
: >- >-• 
c 3 
rux oc. in 

a 4/1 Uj W 

I W fsl W 

0 <£ Ui X 

12 2 UJ 

• X UJ W 

: X o <c 

► ZW w J 

• w CD 

i ul 3 o 
: V) o 02 
3t ^ w 
- 2 inx 

: oo 2 w 
c ur w <2 
1 UJ w X 
: Q. < o uJ 

► in w- o 

► > 2 

1 CkUJ wu. 

2 0 X < 
: w w 3 < 
) X o w in 
X w 

C 111 <C in LL k 

I X 02 o : 

W 2 UJ < 

1 «c a w c 
: in w X : 

c w m X o I 
w w 

1 f ull 

X 

‘ inx I i 

J W < I c 

I O ID O - 
. 2 W X 2 < 
ixm XX 


*00000000000000000000 

JOOOOOOOOOOOOOOOOOOOO 

^000000000000000 00000 
«*ooooooooooooooo 00000 

IvOCOC'OOOOWWW-HWWOOO'CO nOiHUI 
C WWWWWWWWaHWWW 


J—CDOOOOOOOOOOOOOOOOOOOO 

CIUJOOOOOOOOOOOOOOOPOOOO 

llJOOOOOOOOOOOOOOOOOOOOO 

X *^00000000000000000000 

Q owegeg<J^o»oO'n<Hoor-iogi»oo'of^O'tf<g 

Qcs <0 <0 o r- vo «o 0) in ^ ro 

03 S 10 »0 


OOOOOOOOOOOOOOOOOOOO 
I ^ 00 0000 00 OOOOOOOOOOOO 

I 00000000000000000000 

^ co^o:t‘Oor**>oocoinocgowvnoo(7^o 

' wo^iooioinfotnoo'foinioof-awiooino 

I wcgcgminsor^oowcg^ininwwocgm 


D-2 



fO 

lO 


Ui 

o 

a 


in * 

I 

N, ♦ 

xO I 

•H ♦ 

>. I 

r* # 

o 

lU 

H 

< 

O 


* 

I 

* 


10 

a 

UJ 

> 


UI 

o 

o 

X 




>- 

< 



-4 









•> 





r 

K 


O 


o 

o 

o 

r- 


o 

o 

CO 

in 

O' 

ui 



3 


o 

3 

3 

to 

to 

o 

3 

to 



> 

•D 


3 

3) 

3 

in 

J) 


ao 

o 

3 

3 

•H 

O' 


£ 


• 

o 

• 

• 

m 

» 

• 

• 

• 

• 

• 

• 




r- 

+ 

o 

St 

a* 

OJ 

<M 


O) 


O' 



o 

* 

to 

in 

«-l 



to 

to 


r- 

O 


O' 


a. 

o 

o 

w 





Si 



n 

-Vi 

CM 


\n 

• 


h- 











aJ 

s 

4> 


O' 













h- 






3 0-1 


3 

3 

3 



3 


3 

3 3 3 X> 


O' 


0 

1 

X t" 

J) 



II 


3 3 0 


O 

3 

O fO 


to 

3 


3 

0 0 0 3 


r- 


r- 

> 

O 3 

< 





O + • 


O 

m 

Si CO 


ac 

O 


O 

o w in O' 


O' 


3 


2 £ 

.** 





• O' 


* 

» 

» • 


• 

• 


• 

m • • • 


• 


• 


3 in 






3 Si 


o 

it 

4- CM 


CM 



O' 

Si 3 


O' 


3 

-I 

< 

ul 



» 


in O' 


w 


ro 


to 



r- 

ro 


in 


O' 

UJ 

J UI 

X 



CM 


in •• 




» * m 


in 




3- 


CM 


CM - 

a 

u 

H* 

UI 


Si 


o O' 



o 

C3 o 





o 





o 

o 

J < 


3 


II 


O O 3 



3 

3 3 





3 





O 

X 

< X 

O 

H- 


*-• 


•. 3 + 



O » 

i/) » in • 





• 3 

* * * * 




- in 


£ in 

K 

I- 


»- 


O • • CM 


o 

♦ o 

• o •^«>. 



o 


O • 

O O O lO 


to 


N • 

a: 





"y 


O C\J ^ 3 


o w o cr o cr to 


to 

o 


O 3 

O 3 O CM 


c^ 


to 


o z 

H* 

f 




3 to to 


o w m 

Si ?o 


CO 

o 


o m 

o m o o 


■o 


«-4 Q\ 

>- 

z o 

D 

m 


O' 


• if !"• 


• 

• 

• • 


• 

• 


• 

» • • • 


m 


•CM 

•X 


a. 



Z 


a *~t » 3 


o 

3- 

•a- CM 


CM 



(^ 

in -• to 


O 


r- 

-1 

Ui 3 


Uf 


II 


O f'* o 


<-t 


to 


to 



r* 

CM St 


'0 


O' 


J < 

-H 

3 


•»• 


3 • 3 a- 



• 

- - J> 


Si 


* 

» 

■o 

* 

"M 


M - 


t- J 


< 


3 


O + • 



O 

o o 



1^ 

o 

o 





3 

Ij 

r» 

X 

<• 




O CM *-l 



c 

o 3 

to 


to 

o 

o 


O' 


3 

to 

3 

D J 

3 





« 3 3 



o 

j> • in 

CO 


3 •- 

3 

^ o 

» » ^ « 

O' 


•H 

- Si 

;£ 

Z < 

U. 


z 

O' 


O • 30 


3 

• o 

• 3 •r* 

• r-* 

• o 

• o • 

3 3 3 

• 

CM 


•-I • 


m n 





O O CM - 


O^OctOcftOCMtOCMO 


3 lO 

OOOr^tO*^3®l^ 

u 

X 

rii 


2 

II 


3 in »o -0 

O' O *-» 3 

in CO tO :o to o 


O .0 

ooooor-mr'-C' 

«JW 

ui 3 



J3 3* 

• CM 3 3 O 


■ 

• • 

Si 

• 

«n • 


• 

• • • • o 

• 

CM 

• iM 

n 

3 2 

t 



a 

* 

^ *H -O •♦■ 

*!• 

o 

a- 

O- *M 


M 



M 

3 M O' 

M 

3 


•s. 


-t 

•» 


3 

•V 

^ O •‘0 

Cf 


•o 


•o 



O 

Si r- 


0 


O' 


X 






SM • M 




* 

m 

•b 



CM 

-CM 


CM - 

t 

sli 



3 -O d- 

O O 

M 


3 

o o 




3 

o 




3 

Si 

j) 



< -H 


O Si 

XJ 


3 

3 3 

-o 


•o 

3 

3 


O 


O 

J* 

< 

UI 

•-4 


X 

It 


•> O * Cl- 


O 

tn * «n « 

:0 


3 •> 

O 

3 

^ ^ rO 


.O 

- St 


3 

3 


/v» 

j) 

^3 • r** o r» o 

• o 

• o •r«- 

• r'* 

•O 

• 3 ♦ 

O O O O' 

# 

O' 

• 

•H • 


♦ 

3 


3 

M 

t 

3 n 3 o 

3 

3 

-« 3 

t 3 ± o 

*M 

o 

M O 


3 >«• 

3 3 0 0 

3 

0 

o 

M ^ 


•p 


< ♦ 
< I 

u ♦ 


< 

2 


UJ o * o »4 rg 7^ 

j 1-in -H 


• J) • J7 
CM <M 

to to 


3 


1 

# 


r- II - M 

^ 3 - - 


-3 - 3 - 

3 


* 

t 


’ *n o* 

o P-* O O' o 

3 

o r- 


* 


♦ 


P X •* 

3 3 3 3 3 

3 

3 0 0 


1 




X r J- •> 

o t^ ♦ a- •►3 • 

3 - 

3 — .0 — O — 


» 




n - - 3 

• • 0 3 »3 

•3 

• r— • r- • 3 


1 




3^0 

0'^3f^3-aO4'3 

d- O «M to CM 3 


* 




• oJ - - OO' f-r^3-43 

3 

3 to O to O 

3 

1 




CM O 3 • 

3 — • O' • • 

• 

• Si »Si • 


« 




r>.<x--±-4-n-«3 t 

d- 

M M 

J 





X 3 d- O' 

r- 3 3 


■O to 

*3 





t^Ji - II -• 

- O 3 d 

* 

-3 -3 - 






•3 3 3 

3 + • • 3 

3 

O P- f- 

5 





3 jJ - O 

0 0-4-4 3 

3 

0 0-0 

o 





-IMS- 

3-0 - 3 • 

wO - 

3 - 3 - 3 - 






-1 X - M 3 

• 3 O *3 

• 3 

• P» • f-' • 3 






M Z -H -40 

3^ - -3-40d-0 

d 0 M O M 3 

■a 





sI--33M--03-in 

n 

3 0 0 0 0 






tj .< O O • 

3 M O 3 • • 

# 

• 3 •Si • 

< 





i/1 UJ - 

^ d- + + 3 d- 

d- 

-M CM 






< ro .M 

• 3 a 


rO to 

a 





Z X - • 

- d- 0 •- “ - 

* 

- 3 - 3 - 

n 





3 O -O N 

3 r* 3 o 

3 

^ p* r- 

-a 





3 3-d- 

3 CM r- <3 

3 

3 to O 

«» 





II It 3 d- II 

3 II 3 3 II -3 II 3 II 

n II -0 II 3 It 






D n II M 

• 3 3 M ^ 


• D • o • o 

3 



\ 

CM 

2 <X II 

O 3 3 cc ^ < d- uJ d* X M >- -M M 

3 




X »- <J »-• jC 

in • • < :D 

*3 

•3 to 3 0 0 

3 




Si Si S. Z D 

P- O d- 3 “• 

t~4 

-« 3 3 -4 

J 





aS xji X 

^ Ji 

3 

•3 3 3 

3 




^ X r- 





f>J 

00 


o o 

^ • o 

• o • 

O 4> 

o o 


• m • • ^ »:\J *01 

O t) 30 -* -I 
O cO I) 

^ k M • rv) * 

-o 

o VI 

•• o ^ Si 


o 

o 

o 

• o 

3 ' 
O ' 

• 

o 


o 
o 
' o 


O *3 O' • 
o o r- o 
o o--« ^ 

• • • *o 

« ^ ^ 

n sD 


0*0 
o r- 
O' n o o» 

• M • 


. :M 


o 
o o 
O 4> 
:3 
• II 
0 > 


r 


< 


0 
s 

^ *. •►O' 

o o%\) • 

3 3 -M 3s 

3 3 O' n 

» • • CM 

»0 O O' 

^ -o o D O' 

•- \J -M 

O- 


J) 

II ri II tn 
^ • 
u> — J -o 
uj.n 

JD 3 


^ 'O 
3* 

n ---o 

o • ^ • 

3 n M r- 

0 O 3 > 

• g • CM 
O /) •> 

o 

- 3 

(\J O 
O' • 

H 3 

r^ 

O' ^ 3 

N z 

II zi 

r ^ 


-o 
II II 

m Z 

J 3 X I 
•iJ J) S ( 
3CM J 


D-3 



it 

m 


ui 

o 


«n 

r* 

V 

vO 

»-l 

f- 

O 

UJ 

o 


» 

* 

I 

♦ 

I 

* 


llJ 


X 

ui 

> 


III 

wj 


in 

ui 

o 

o 

X 

ct 

LJ 

. >« 

< 


3 


u 

u. 

in 

X 

\ 

c 

m 

z 


o> 

•H 

cc 

>- 

in 


txl 

o 

X 

3 

o 

<n 

:a 

z 

< 

UI 

sn 


Q 

::> 

o 

o 


X w 

o r> 

z X 
zj m 


» 

t 

* 


z 

o 

w 

cn 

X 

UJ 

> 

Ui 

a 

o 

X 

X 

UJ 

>- 

< 

-J 


-J 

D 

X 

o 

u. 

\n 

X 

N» 

«« 

UI 

z 


in X 

X 
Ui o 

u z 
< 

X 

in 


Ui M 


<VJ 

X 

z 


X ♦ 

X I 

c ♦ 

X I 
Ul » 

t 

« 

« 

I 

» 

¥ 

t 

* 


'xl 

X 

< 

2 


0 


0 

0 

0 

r- 

r- 

0 

0 

0 


0 

0 

0 

rt 

rt 

0 

0 

0 

O' 

0 

if) 

in 

CO 

00 

0 

0 

• 

0 

• 

• 

• 

■ 

• 

• 

— 

f- 

♦ 

0 

4 

4 

cy 

cy 


CO 

rt 

rt 

-4 



rt 

rt 


r* 

31 

O' 




m 

m 




II 

CC 

fNi 


tn 

r- 

00 CO 

O If) 

+ • 
o ^ 
OJ 
cy 


o 

o 

o 


o 

u> 


o 

in 


fO 

CO 

• 

lO 


r- o 

KJ- O 
00 O 

» • 

cy 

fO 


o 

o 

o 

• 

O' 

r- 


U Ul 

X 


-4 

m — r- - 

- - -in 

in 




— 

0 

w 

Ul 

^o 

0 CO O' 

000 





0 

_j < 



11 

0 33 0 

000 





0 

<£ O: 

0 

w 


- 0 0 

- 0 -in - in - 

— 




- 0 

X in 

1 - 

H- 

>- 

0 • »o 

0 *o »o »r- 

r— 


0 


0 • 

QC 


3 

z 

0 cy rt 4 

0 w 0 4 0 4 rt 

rt 


0 


0 3 ) 

0 X 

f- 

X 


0 rt r- 

0 -4 in in <0 

CO 


0 


0 ul 

z 0 

3 

in 

O' 

•4 CO 

• • • • 

• 


m 


• 

z 

a 


33 

0 - 4 - 3 ) 

0 4 4 CM 

cy 




O' 

Ul 3 

z 

UJ 

II 

0 r- r- 

-4 rt 

rt 





--J < 

w 

0 

w 

in -0 4 

- - -33 

-in 

— 


— 

— 

I- ^ 


< 

U> 

0 + * 

000 

r- 

r- 


0 

0 

J- 

IT 

X 

z 

om ry 

000 

TO 

10 


0 

0 

3 J 

0 

in 


- 0 GO 

-0 -in -in - 

CO - 

CO 


0 

- 0 

X c 

U. 


xc? 

0 »C 0 

0 -o *0 -P- 

»r- 

• 0 

• 

0 • 


0^0 oom •» •-o o ; 0 - o ^ fO oj lo oi o 
zii «^otninooa'o*Hvn m <0 to :o o 
Z 5 i-<:r *cyiooo » • • *m •*n 


•OOJt'-H 

c\i «^c\f o in *> 


-Jrt 5 *- 

< r-l 

X tl 

X in 

O NI :f 

z z 

UIO - 
-I 

»- M - 
H* Ul ^ 
D X 
X Z :» 

Ul ^ 

O ^ 
• Ul -> 

Oi U O ^ 

< * • 

tl^O :t 

i-<n n 
o 0:3 

O UJ o 
^CM s: 
W CL -^rsi 
rj Xrtl 

o X O fO 
m Ul •• 

^ rt CM 
X r • 

o aj m r- 

sO \D ^ 
II II o ::r 
OU> II 
CMZ <(X II 
£ w O 1-4 X 
«t UIX 'ICO 
z tti < m < 

m i- z «-•)- 


o 

o 


CM O' 


^ O » CO 
O •h* w 
O i> O JO 
o ro 4- • 

• ^ O -fl 

o w ni 
m O' 

OJ •.:S' • 

or- ® 

O :3“ O 

o ■-< + 

o « • 

o co 'O r- 

O O' ^ 

• O CO 

ct — * » m 

o* r- in 

-I • o r- 

O 'I' • 
o -o -M 

o *.0 

O O CM 

o o r» tJ 
»oc0 o 
^ T>» ^ 

»tn 

nO 

O <M 
O CO 
ti o II <y 
fSl *3 O 
o in 
m • 
r- ^ 


O' "* o 

•HO • 

O 

« O «-l 

in • 

O' 


O' o 
o o 

4- •• O 

in o • 

fO O *4 

m o 

lO • 

O' ® 
o 

• o 

m o 
•*0 

o • 

•. o w 
0 0-4 

O • 

+ >0 
«n 

in a 
KJ o 
ON II O 

r- z • 
o X 

1^ 01 


^ ir CM 
K> 

» •• m 

o o 
o o 
—\xy * to - 
o • o »r- 
o 4- o jJ“ *0 
j) in CO 
• • • 
if o- <y 

K1 

* 

o o 
o o 
-*n •►•JJ •- 
o • o •r* 

O if O 4“ >0 

^ in 33 


•n 

*in 

r- r- 
fO m 
CO •. CO 

• r^ • 

M -o CM 
10 CO fO 

► m *10 

CM 

lO 

• tn •- 
r- r- 
lO -o 
CO *> CO 

• f- • 

cy.-n cvj 
fO 03 /o 


in • J 3 
4 * 4 * cy «y 

fO *0 

1. •►m -in - 

o o 
o o 
-•in -lO - 
o *o -r- 
O if o 4' fO 
n J 7 JO 


f'- r- 
m Ki 
> (O • CO 


if ^ 04 

lO 

- - in 

o o 
o o 
II in II in It 
iC * aC • O 
< if Ui 4* X 

000 

w *-• 

in u) <n 


OJ -V) CM 

K> JO »0 

m •in 
CM 
to 

- n — 

r- r- 
fO to 
m It CO 

• o • 

01 >- rvi 
lO o»o 
in Win 

in 


CO tn O' 

.•o f- 

O W O' 


*H 0> 

o m 
•O CM 


o o o »n 
0 0 0 33 
0-4 O O' 


O O O fO 

o O O CM 

0000 


33 

O' 

OJ 


33 

r- 

CO 


MO O' 33 

m m O' 

4 - t\j cvi - 

o 
CO 

• in 


to 

o> 

ro 


O to 
O 3) 


<y 

3 ) 


m -4 in 10 o r- 

CM -4 if 33 o» 

lO - CM - ai - 

r- r- 33 

O' 33 K 1 

• •.*. »o> -w - m 

OOOf- -CM --4 • 

ooor-fO«-»ODor- 
oooocor^mr-O' 

• • • «o *01 * cy 

o CM o O' ai o r- 

m rtf- 33 o 


o 
o o 

-o -o 

3 * 0 * 

o o r- 

O O 33 

<y 

CO 

•- 

> o 
o o 

-o - o 
o -o • 

O O CO 

O O 33 


cy - 

cy - 

cy - 

4 

0 

in 

rt 

3) 

4 

-rt 

-rt 

- 4 


0000 • CT* • *4 • 
oooooortvr- 
oooifwwininO' 
• •• —•O 
o 4 * o to »4 -4 r- 

O 33 CO sO O' 

.CM -CM 


O 
o o 
- 0-0 

o -o • 

O 0-4 

o o r- 


o o 
o o 
o o 

o o 

0 -4 

01 


o 

CO 

o o 
o o 
It O IJ o 

o • 

IM <C 3> 
H- r- 
UJ 


o o 
o o 
o o 
— • 
o o in 
»o rt 

O 4 
vO 

to 

M II 

II W IH 

>00 
-«c 


• f- 

to 

- -r- 

O O' • 
o r- o 
o -4 4 

• *rt 

o r- w 
o 33 
-4 - 

•D 
4- 
-.O' 
O CM • 
O CM O' 
O O' 33 

* • ^ 
0 0-4 
33 rt 


tn 

tt II in 
-4 X • 
H -Jrt 
ujtn 
d o 


r- rt 
n cy 

- o - in 
33 • n • 
o 33 r- f- , 

O' in o O' 

»cy *c\i 
o r- 

33 O' 

<M - cy - 
4 rt 
4 O' 

— 33 -1 rt 
o • 4 • 

o 33 <y r- 
33 33 3) O' 

• 'M *01 
o m - 
n O' o 
<y •» <y • O 
4 CM O 
rt O' • 

It r- II -4 o 
>* « <^ -O' 

-j 33 X p-r- 
UJU3 X O' W Q 

o ai *U <M z 
H- If ul 

X ^ 


D-4 


CLOUD Rise AND SOURCE DISTRIBUTION NASA/MSFC MULTILAYER MODELS V -5 DATE 07 / 16/75 


in 

fO 


X 

1- 



o 


o 

o 


o 


r- 

P- 

o 

o * 



CO 

in 

O' 

UJ 

-1 



o 


o 

o 


o 


to 

to 

o 

o 



fO 

r- 

r- 

> 

3 



o 

O' 

o 

m 


in 


00 

CO 

o 

o 



o 


O' 


£ 



• 

3 

* 

* 


• 


m 

• 

• 

• 



• 

• 

• 





r- 

+ 

o 



it 


CM 

CM 


CO 



«4 

O' 

M9 


O 



fO 

o 

-4 





fO 

to 


r- 



o 

in 

O' 


u. 


o 

S} 

f- 






in 

in 





xO 

CM 

CM 


in 


• 


O' 















UJ 

z 


CO 


00 















_f 

•s 


«-4 


- sO 






» 

* 








1- 

< 



O 

CO O 

o 

o 


o 


r- 

r- 

o 

o 

o 

O 

o in 

c^ 

vO 

X H> 

m 


II 

o 

o o 

o 

o 


o 


lO 

to 

o 

o 

o 

•H 

o ^s> 

r- 

r- 

u => 

< 



o 

+ • 

o 

tn 


in 


<0 

00 

o 

o 

in 

O 

O O' 

O' 

CO 

Z X 

z 


o: 

• 

CO CM 

# 

• 


9- 


• 

• 

• 

* 

• 

• 

• * 

• 

• 

o tn 



M 

CO 

w 

o 



it 


CM 

CM 


O' 



*4 CO 

O' 

^o 


UJ 



in 

o 

*4 





lO 

to 


r- 



to 

in 

O' 

-1 ut 

X 



tn 

* vO — 





— 

tn 

in 


— 




CM 

CM • 

L> 

K 

Ul 

vO 


O CD O' 

o 


o 


o 




o 





O 

J < 


-1 

II 


0x0 0 

o 


o 


o 




o 





00 

< Q> 

O 

►- 

•-4 


O CO + 

— o 


in 

* 

in 

— 

— 


• o 

* 

* 

— — 

* 

• in 

£ V) 

H 

►- 

H- 

o 

• *00 

o *o 

• o 

•f- 

N 

o 

o • 

o 

o 

OlO 

lO 

r- • 

ttl 


3 

Z 

o CM ro rO 

o w o d* o ^ to 

to 

o 

O M> 

o 

lO 

o CM 

O' 

to t- 

O X 

f- 

X 


o to O' 

o-^ m 


in 


<o 

CO 

o 

OlO 

in 

o 

o o 

ro 

-rH O' 

Z o 

D 

in 

O' 

• dr M) 

• 

• 


• 


• 

» 

• 

• 

• 

* 

• • 

• 

• CM 

z 

OL 


(O 

O <-4 - CM 

o 

d- 


it 


CM 

CM 


O' 

CM 


in rO 

o 

P— 

U1 3 

z 

Ul 

It 

o 

f- lO 

*4 





to 

to 


p- 



if 

xO 

C^ 

-* < 

»-* 

v-> 

*-• 

lO 

- O *4 

— 





in 

- J> 

* 

» — 



to 

- CM 

■w 

- CM - 

r>^ rn 


o + • o 

O (O fO o 

» o >o o 

X O' O •h- o • < 

U\0 O O Cf o ( 

z II »o m o (X) o«-i L 

*(> 4(000 *■ 

< a *-i 3 + + o 
^ Z :t<0 

«liO ;f o o o 

* o oo o- 

XII 4 *«- 0 *>(Tvd* -“O 
cc (/) *o *h- m in o •< 
O*' 4 :t 3 J^O' 03 O*H 
Z Z •.O'O -f • 


O 1 ^ o o 

o ro Ki o o 

kin » CO ^<0 

> *h“ *o *o • 


r- (O 

O' vO lO 
» -nf) 

OOOt^ *fU **-t • 


O . 

it to 

CM 

to 

CM 

o 

o 

to 

o 

o 

o 

P- 

to 


o 

CO 

p* 

in 

00 

fO 

o> 

to 

o 

o 

vO 

o 


o 

o 

CO 

f- 

in 

p- 

O' 

• 

* 

in 

• 

m 

* 

• 


• 

• 

• 

• 

o 

• 

CM 

• 

CM 

it 

CM 


CM 


CM 


Si 


o 

O' 


o 


1 - 



to 


to 



:0 




•-I 

r- 


■o 


O' 



-in 

— 

tn 

— 


* 





CM 

— 

CM 

* 

CM 



II 

in 

o 0- 

X z 
■j} * 
o 

* m * 

c \4 U O 
< *. 
0. o 
>-in •* 
o o 
o jJ - 

-104 

•^a. •> 
CM r :t 
< * 
0 X 0 
inuj 
X lO 
XX - 
O O) lO 

^ ^ «k 

II II o 
o«n II 
<M Z <0. 
Z H> O w 


iii o •> o ^ vO O' 

-I r* 

I- II -04 - 3 ) * - 
m (I* o o GO <> 

3 0 . - o JO oo 
XZjJ'-O 0 »+-** •• 
lO - *0 • • j±- O' o 
OjfOC0»0»0OO 
m * * o O' »o iTi o 
uojf^o inoj* 
<( *. -inoco 

OL O ;J- O' p- cO o 
U) -11^ -od-iO 


> o r- 

> O lO K> 

> — in — co - eo ^ 

• o •r- ♦o 

• O ;t *0 M -o M o 

«n CO to <0 lO o 

• *m *10 • 

^ CM (VJ 


o o ^ o m 

o O K) vO ^ 

0—0 ^ — — fO— lO— ^ 

• O • OOOO' •O' • *-• • 

o h- oooGooonMP- 

o M) o dP o ^ in O' 

• ••••X)*CM*CM 


CM 

CM 

CM 

o 

in CO *4 *4 

p- 

to 

to 

CO 

*4 

CM CO M> 

O' 

■s% 

- m • 

* 


- CM 

- CM 


► o r- f- o o 

> O fO (O O O 

) — lO — ^ 0 — o — o 

• o *r- —o »o • 
■ O KJ OJ lO <MO o a> 

in ca^ootoo o^o 

• «vn •m • • 

:J- <M CM -• 

lO lO GO 

k — — 


f- f- K> 

to J> CM 

— — — — f— — O — iO 

OOOO* •o *%o • 

ooor*ooj>p-r- 

OOO^^O'O'vrtOO' 

• • **10 »CM *CM 

o o P- -« O P' 

CM o >0 sO O' 

*-• - CM - CM - 


< inz ^ 

Z ul<in < 

w *- z 


O o + • • o 
o o ^ »o .o o 
X - o -O 

MO *lO o • 

►-40 O -I - - o «-< 

-OO0'3330O*H 
^ • O O O • 

(O ^ 4 - M> 

CM *p-o 
* — O' O GO — 

P*. O ^ o 3 

O f-co o 

^ II o II ofo II o 
fM • 3 o ;0 3 ^ • 

II O 0 ' 0 'X ^4 
5 ^ lO • • <*^ 


O O 

o 

P* f- o 

o 

xO 

d* 

lO 

o o 

o 

0-0 0 

o 

d* 

d 

c^ 

*»o -in 

- J> 

— ^ — > dJ — o 

- o 

•• «k - •iO' 

- xO 

- o 


► o *o •r— •p- *o *o • 

IO^OJf»0<M0CMO o-^ 

n ^-oujjoo or- 

• ♦ •in»m* • 

^ ^ CM CM O ( 

K) to d) 

I — •m—m— — — c 

► o o r«- f- o o c 

» o o po to o O V, 

t II in If in II CO II CO II o II o 

• *:^ *0 *0 »0 *a^ ♦! 

► < 4 'U>^XCM>-CMIM <M>; 

lO O O to 3 iO O H-r- « 

Hv uj 2 

u) lo «n %n i/i X " 


) O O CM »o • • 

«3 3 <MO'OinMf- 
>ooo'j)LOin OS' 

» • *»CM *CM •CM 

I CM O O' «H O Si — 

in fO no O' o 

- rM - CM - 3 

^ it CM o 

m to O' • 

I II II II »n II r- II 3 

• ')C •O' 

ko»--»»o_j'OQ-r-p- 
ujin ujin X O' ^ I 

O O O CM ul CM 


D-5 


>- ^ 

-JO 


yO 

to 


U1 

O 

< 

a 


in 

r- 

X 

yO 

V 

o 

Li 

K 

< 

O 


«n 

I 

> 


•o 

»o 

z: 

a: 

ui 

> 


Z) 

z 


in 

X 

>s. 

<C 

I/I 

< 

2 

2 

O 


z> 

aa 


I/) 


UJ 

o 

a: 

3 

O 

m 

3 

Z 


111 

UY 


a 

3 

O 

-I 

u 


I 

* 

I 

« 

a 

« 



o o 

3 3 3 3 



3 33 3 



3 GOO' 3 
O' *4d“ O 


X X 

• • • • 


J UI 

O' 3 O' 

o 

UI >- 

fO 3 3 

o 

O til 

3 

in 

X 


• 



if 


o o o o 
o o o o 


O X 

o o o o 


N Ui 

o o o o 

in 

o >- 

• • • • 

•-I 

»H UI 


in 

£ 



in o»^r^rocn wr^tnoio*-iGOoc\joj^m o 
o o^oo0'C7'v0ininr0o<0f^i0*-*inv0K) 


I- 

UJ 

X 


O X 

>- u 


oj fgcsjc\jcMc\jc\ 4 CN 4 <\joj(\i<\ 4 cuc\j<Mcgavr*-w 

fOroro»oioioiOfOio»oK)fOK)>or<)cofOh-in 
CO (ococo<x)co<ooo(oco(Dcoa><ococoinr^ro 


» 

I 

* 

I 

* 

I 

« 


til 

-1 

u 


X 

Ui 

> 


< CL 
X tn 
tr 

o X 


< * 
X I 
UI* 

« 

» 

I 

» 

I 

« 

I 

» 



o 

»-i til 

3 

3 3 

3 

3 3 

3 3 3 

3 

3 

3 

3 3 

3 

3 


-H 

r* 


• 

in X 

iO 

310 

3 

3 3 

33 3 

3 

3 

3 

3 3 

3 

3 

O 

3 

m 


3 


m 

m m 

m 

mm 

mmm 

m 

m 

m 

mm 

Ip 

m 

m 

if 

3 

O' 
















if 

II 


3 

3 3 

3 

33 

333 

3 

3 

3 

3 3 

3 

3 

O' 

r- 

•4 

•4 

X 



l^f“ 

r- 

r- 


f- 


r** 

r-f- 


r- 

if 

3 

3 

o 

X 

O X 

fO 

to to 

3 

3 3 

3 33 

3 

3 

3 

3 3 

3 

3 

3 

f* 

m 

o 

ISI 

X UI 

3 

3 3 

3 

3 3 

33 3 

3 

3 

3 

GO 3 

3 

3 

m 

1^ 

3 

• 


o >— 

• 

• • 

• 

• • 

• • • 

• 

» 

• 

• * 

• 

• 

* 

• 

• 


H- 

»-» Ui 

3 

33 

3 

3 3 

3'3 3 

3 

3 

3 

3 3 

3 

3 

r** 




X 

in X 

lO 

3 to 

3 

33 

3 3 3 

3 

3 

3 

3 3 

3 

3 

o 

3 

m 


o 

►H 


m 

mm 

m 

mm 

mmm 

m 

m 

m 

mm 

m 

m 

m 

d- 

3 


UI 

X o 


o 

o o 

o 

o o 

o oo 

o 

o 

o 

o o 

o 

o 

o 

o 

o 

in 

o 

X 

o 

o o 

o 

oo 

o o O 

o 

o 

o 

o o 

o 

o 

o 

o 

o 


»- m 

UJ o 

o 

o o 

o 

o o 

o o o 

o 

o 

o 

o o 

o 

o 

O 

o 

o 


2 • 

O UI 

in 

m m 

m 

mm 

mmm 

m 

m 

m 

•H *4 

— 1 

•4 

•s 

-H 



« 

I 

♦ 


if 


-I 

UJ 

o 

o 


X 

Ui 

UJ 

X 

V 

X 

in 

Ui 

UJ 

X 


UJ if 
X 
UJ 
O X 
O 3 

o in tn 

• < w 

UI 

X X 
X 
UJNI 
X 


»-• Q 

in ^ 


X ^ 

< o 

O Ui 

•-» a 
in 





o 

in f-v- 


ui 

X 


UI 

< 



o 


Q 

»- 


1- 

X 


3 

m < *oi 

3 

X >- 



r^2 

X -1 

o a 

in 


3 w 

UJ UJ o 

O 

2 X 

X 


• 

1- -J 2 

3 

a m 

UI 


n- X 

UI o < 

-J 

< 

1- 


if X 

X 2 

<C 

-1 UI 

Ul 

o 

•d- UI 

<c 2 


u 

X 

o 


2 O 


-J < 

< 

o 


►-< X 



X 

c 

X 


o „ 
o in < ^ 3 c 

o »-MM cn CM X > 
X UJ 


UJ X *-• X 
X OQ O X > > 
►13 <X O I- UI UI 
O UiOX UI X X oo 
3 X-J 0>- O 
O U C X r-i 3 3 
33 33 UD 3 XX 


^f-f-if^sf^ififif^ifif 


O 33000000000000003 3 
O OOOOOOOOOOOOOOOOOQ 
O 000000000000000030 

tn tninininininininin;Oin*-i^^3*^^ 


r- 

o o o 
+ + + 
in -« CO 
^ vO n 
in 1 ^ 

•H CO ^ 
O' oo 
;f CO CO 
CO C7' O 


00 O CO 
0 0 3 
+ + + 
3 r* d* 
W OfO 
0^10 
o r-m 
3 otn 

0 0 3 
3 Oh ^ 


3 Oh O' 
O OO 
f + + 
^30 

O Ch O' 
O 3 3 

a>c\i*o 

into o 

10^0 


3 Ch O' 
0 3 0 
•f + + 
O O' 
3 O if 

0 m r^ 
3 CM o 
.3 O if 

mom 

01 fO 3 


O' O' or* 

0 0 3 
+ + + 
O O' 3 
3mm 
03 O' 
c^m CM 

3 CM w 

m o O' 

or- -I 


O' o o o 
o <H 

+ + f + 
K> m h- CM 
rO 3 3 h- 

if ^ m 

O if CM O' 
3 >0 CM O' 
«H ^ O' O 

3 f- in ^ 


0>K)»OK)^Oh*fO<M3K) 3 K> K) r- ^ -I ^ «H 


2 U 

o 

•H O U2 UI3 J 




2 

UI 

0'uiuiifin>><ui 

1 


o 

UI a 

)- 

r- m >- i-i < 



o 

-1 < 

< 

m -4 o J J a o 

3 


o 

f- -1 

-1 

2 X UJ if 2 2 


3 

o 

w 

a 

Q X »-c 

X 

O 

o 

a -J 

o 

m m >-0-1033: 

f- 

O 

o 

X < 

-1 

X »-i2 2 X UJ 

o 


o 

in X 

< 

UJ o UJ a X UJ UJ 

2 


• 

X 

o 

1- m *-i »-i X o o X X 

3 


o 

UI o 


Ui XI- 31- X l-l- l- 

X 



O 2 

» 

X UI <X»-I »- 

1— 



< 

1 

1-3 < 3 X O X X 

m 



X 


2 X»-« 1-30 0 



o 

m 


w X -IX UI *-« 

'aJ 


4 



»-»o in 3 UJ 2 2 

o 


3 

Ul 


1-2 3 a xo o 

X 


d* 

_j 


X »-• < 1— UJ f- M k-» 

a 


3 

X » 


5 1- 2 X n 1- 1- 

o 


3 

X 1 


»-i>-m3oai-<< 

J) 

o 



O 000000300000000000 


d- 


r* 3 3 3 
o o o o 
+ + f f 
m o m f- 

3 CM 3 r- 
3 CM CM if 
m r«* 3 3 
m 3 CM m 
to m o m 
if o m 


3 O' O' 3 O' 
OO o o o 
+ •♦■+ + + 
3 o to 3 m 
o d* O' in lO 
cMf^ m m m 
f 3 *o o 
3 3 3 O' C^ 
^ f *3 O 

3 d* m o 3 


O' O' O' 
o o o 
+ f f 
d- d* CM 

mom 
O' cMd- 
r- if 
O' cMr«. 
^ OCM 
d* d* O' 


O' O' o O O' 
o o-« ^ o 
4- f + f f 
CM 3 3 m 3 
o O' 3 o 
CM if ^ ro 
3 f I'- O' O' 
m if 3 m 
3 «o m m 

3 d* K) M f* 


r-KjcMfO ^cMmfM- 4 f-<MmcMCMin-H^-« 3 


oi- o X ^ a a 


4 

O-f *-#0-l>-22 


o 

3 0 31-)- X < < 


3 

I- X X)-X «C)- >- 


3 

m j»-*3ooomm 


r- 

a <1- .> 3 »- 2 

-J 

3 

0 3 ail 

o 

3 

a X 1 1 f 1 o I 

X 

if 


f- f- r^ 
o oo 
+ f 


< 3.x X 

I M I- X <C aim 
t aOOOI33X 
X<X3I- *-#«-txI 

UJ X 
X UJ X 
^ rsl >- O 
X <H- 
<1 -I 
X 

X X 
UI UI • 

>- >- o 
< < 2 


3 r ^3 
o o o 
+ + f 
n- m 3 
in 3 3 
O' f*» !>• 
*4 CM 
r- m m 
K> 3 
m m 


3 0'3 
o o o 

+ + f 
3 cMm 
-4 3 ^ 
3 »or^ 
m I'- O' 
if 3 -4 
oo W 
#o if o 


3 O' O' 
O O O 
+ + f 
3 f^if 
h- r- CM 
o f- CM 
O' 3 
o ^ ^- 
cM 3 r- 
o if o 


O' C7' O' 
o o o 
f + f 
CO ro n 
o r** ^ 
CM ro m 
3 3 CM 
d- CM O' 
O 3 CM 
<>0 3 


3 O' O' O' 
o o o o 
+ + 4 + 
C^ O' O O 
fO m O' <) 
3 o ^4 3 
O' «-i 3 
K> d- O d- 
3 O' m 3 
3 *4 O' 


^ ^«Hw 3 ^io«^ 0 'd‘WK).^w*o 3 f*r^d‘ 

o oooooooooooooooooo 
ooooooooooooooooooo 
oooooooooooooooooo o 

ifOifo 3 r^oo 3 mo 3 or*- 3 ooo'o 
O'in 3 oinfomoo*mmf 0 of'^*-iinomo 
«Hi 33 mm 3 r^oo^ 3 ifininr^r^o 3 'n 
W«- 4 W^<H» 4 ' 4 W*-I 333 

W 3 ' 0 d‘ 33 r^ 30 'O* 43 ' 0 d*tn 3 N 30 ' 


D-6 


CUOUO Kise ANO SOURCt OISTRiaUTlOR NASA/MSFC MULTILAYER MODELS V-5 DATE 07/16/75 PAGE 

example SPACE SHUTTLE NORMAL LAUNCH *-*.*-♦ 

NORMAL LAUNCH SPACE SHUTTLE VEHICLE *-*-•-* 

namelist NAM2 for INPUT TO THE NASA/MSFC MULTILAYER MODEL VERSION 5 


IQ 


CO 



^ « 


o o o 

tl 

o o o 


m 

ir 

• • • 

N 

in 





txJ 


J II 


1- w 

% * ^ 

h- 1- 

o o o 

3 2 

o o o 

X » 

o in o 

1/) O' 

• • • 

«o 

in w 

llj II 

rsj 

u *-» 


<c o 


a 2 


lO 

^ 4k ^ 

XO' 

o o o 

O vO 

o o o 

2 II 

o o o 

3 >-• 

• • • 

< 3 

o ^ :g 

-1 2 

in 



-1 rg 


< 


X II 

* A » 

q: (O 

o o o 

O M 

r> 3 o 

2 2 

o o o 




CO 

#o 


lU 






o 






< 






CL 






in 

« 






1 





N. 

# 





vO 

1 






« 





N. 

1 






ft 





o 








» 




ui 


« 












in 




o 


2 






O 











« 

to 





1 

£ 





♦ 

UI 





1 

> 





♦ 






1 

u 





« 

UJ 






o 






o 






£ 






£ 






UI 





111 

>- 





-1 

< 





u 

J 





—• 






X 

>- 





UI 

J 





>- 







£ 






U 

• 





U. 

.o 





to 

• 




UJ 

s 

CO 





V. 



^ ^ ^ 


1- 

< 



o o o 

1 

X 

to 

11 


o o o 

> 

u o 

< 

yc 


O r-» o 


2 X 

2 

£ 


• • • 

10 

Z> to 


Nl 


U) -1 


< 

UI 




111 

Ji UJ 

X 

•-I 



Q 

u 

f- 

Ui to 



O 

-J < 


-J II 



£ 

< CL 

o 

J- •-» 


« 


£ UV 

K 

H- I- 


o o o 

q: 

£ 


3 2 


o o o 

u 

O X 

J~ 

X 


o o o 

>- 

2 o 

3 

to O' 


• • • 

< 

2 

a. 

sO 


in ^ in 


Ui => 

2 

UJ II 


CM w 

rH 

-4 < 


u ►-» 



H* 

»- -J 


< u 



-1 


£ 

CL 2 




3 -J 

O 

tO 


^ ^ * 

£ 

X <c 

U. 

XO' 


o o o 


to £ 


CJ vO 


o o o 

o 

£ 


2 II 


o o o 

u. 

UJ o 

£ 

Z>p-t 


• * • 

n 

O 2 

< 

< 2 


O CM o 

£ 

< 

2 

-I 2 


n to 


cu 






to 

H- 

JC\J 



to 


to 

< 



< 

UJ 


£ 11 


« » » 

2 


-J 

£ to 


o o o 


£ * 

.iJ 

O M 


0 3 0 

Z 

£ 1 

£ 

2 2 


o o o 

o 

< ♦ 




• • • 

•H 

X 1 


UI o 


0^0 

h- 

UJ * 


-1 


o vO 

3 

1 

♦ 

»- 11 



2Q 

# 

I 

to 



►H 

* 

♦ 

Z>tL 



K 

1 


£2 - ^ ^ - 


« ^ ^ * 

^- 

» 


J 1 o o o o o 


o o o ^o 

10 

1 


o o o o o o 


o o o r- 


» 


• Ul«> o OLO*n o 


o o o o 

a 

1 


vMOO • ••• • 


• • • • 


» 


>0 


o o o r* 

Ui 



a. o fO fO o -1 in 


0-4 0 O' 

u 



►-to • 11 ^ 


CM CM 

£ 



L> O O ^ + 



D 



OyJ *-0 *0 



O 



JM£ 


^ ^ * 

o o o O' in ^ 




— 1 0- •■Mo^ooor-r^ioo 





MXoj^or^oo o»o»ono 


o o o M n CM 




< •>. *»»Od'OJ>i)33T5-40 


O O O CM ^ 4) 




O X O rO Td • • •• • •• 



< 



touJ r*» ^\0<^^cM<vjK>oooj>o*0XMn •> 




aC lOvM ♦ lOOcOd) 

• oi*>ino»oo'o 

ui 



xx«-* -< in If) ^ 

O ^ ^ OJ ^ o 

lO 



O CO ro r- 

'O 

' 4* CM 




*0 sD 

K> 

» 

£ 



11 11 O ^ II It II II II It II II II 


II It II II tl fl CO 




OtOll NJ<35^jC^OOOj^II'-«-»^-«X>-X»0 

O 



C\I2<Q.II £<UX>-M<>*Ca(Jh->-l-IQ-0 

ro 



:e>-o>^x <ooooo)-< 

U1UJ£ Q 

o 



<tO££:3 aa oj £ 

aoul , 2 




2Ul<tO< OtOi/)tO\OtOX*-» 

II UI 

u> 



tA^^Z ^ H- h- 


X lO 


D-8 



CuOUD rise and source DISTRIBUTION NASA/MbFC MULTIUATER MODELS V-5 DATE 07/16/75 PAGE 


to 


* 

I 

# 

1 

* 

I 



» 





1 





» 





m 





z 





0 





»-< 




* 

to 




1 

a: 




# 

UJ 




1 

> 




* 





1 





* 

Ul 





a 





0 





X 





a: 





Ul 




lij 

>- 




-1 

< 




0 

-J 




r-f 





X 





Ul 

-1 




> 

3 





X 





0 





U. 

0 




to 

m 



Ul 

X 

(O 



-J 

N. 



^ ^ 

y- 




000 

X *- 

to 

II 


0-4 0 

0 3 

< 

X 


moo 

Z X 

z 

q: 


• • • 

3 to 


M 


—1 

< 

tu 




•I Uf 

X 

•H 



0 

1 - 

Ul m 



J < 


j II 



< £L 

0 

1— >-• 


^ ^ 

X m 

h- 

1- H- 


000 

cr 


3 Z 


0 m 0 

0 X 

*- 

X 


moo 

z u 

3 

in 0^ 


• • • 

z 

a 

sO 


CM tn 

Ul 3 

z 

Ul II 



J < 


U *-t 



1 - -J 


< u 



K 

Qi 

CL Z 



3 J 

0 

to 


* ^ 

X < 

u. 

X 01 


000 

in X 


0 %0 


000 

QC 

(M 

Z II 


0-4 0 

Ul 0 

X 

3 M 


• « • 

0 z 

< 

< 3 


m 0 

< 

z 

-1 Z 


-4 

CL 





to 

1 - 

-J CM 




in 

< 



Ul 

1-4 

X II 


^ 

_J 

-J 

cr in 


000 

a. » 

Jj 

0 -M 


000 

£ 1 

X 

Z Z 


0 4 0 

< ♦ 

< 



• • • 

X 1 


UJ 0 


0 m 

Ul * 


-1 


W CM 

1 

» 

1 - 11 



« 

1 

1- in 



» 

* 

3 X 



1 


X z - • • • 


» » 

» 


m •• 0 000 

0 

000 

1 


0 0 000 

0 

000 

» 


• Ul ^ 0 0 m m 

0 

000 

1 


(MOO ♦ • » • 

• 

• • • 

# 


r^< m.'M •.^ 44 ’ 

sO 

0-4 0 



0. 0 10 fO 0 -H 

m 

CM 0 



»— m •> II 4 - 
0 0 + 

0 Ul 0 ^ 


-4 



^ 0 . •* Ni 0 r- 0 0 0 lO 

0 

0 0 0 ( 



CM 5: 4 * ^0 O'OOOfO'Ono 

0 0 0 ( 



<c ^ ^o-Homm^) 3 D^o 

0 0 0 ( 


o 

r- 

O' 


a X o ro 
to uj •> r* 

ill n <\j 

X X 


O CO 

-o 
II It 
O l/> 
OJ 2: < 
X I- u 
< (/) T 
i: uj «« 

I- 2: 


to r*. 

it 

o ZJ- II 
II rsl 
a. It 

i-« ^ 

3 

lO < 


^O^^t<\IC^J»000 
lo m a) cQ • 

m ^ o 
o 
^£} 

II II II II II II II 
X O O O II 
a: < uj X >- M o 

3t0v0t0vAt0 X *-* 


o oi o r<i min •> 

4“ m O fO £^ o 
^ fSI C4 o 
4 - oi 

II If II II II ft CO 
i-i X >- --^ m 

-JOlo 
lU lii X W 3 
30Ui Z 
►- II *jJ 

X tfl 


D-9 
















-JO 

» 




UJUl 

Ul 




O 3 

lO 





CM 

o 





it 




lO 





CM 




xcc 

• 




•lUl 

fO 




o lUK 

o 

LJ 



O OUJ 

w 

O 



U) X 


< 



• *«# 


ql 









o 









0(C 

o 




1st Ul 

w 




irt OK 

• 




»-• !-• UJ 

fO 

lO 

» 


lAX 

CO 


1 


(/) 

w 

\ 

* 


cc 



1 


ui 



« 


1- 

CM 


1 


UJ ^ 

K 


* 


X otc 

to 

o 



KUl 

CO 




o OK 

• 

hi 



O W4 UJ 

CM 

!• 



• lOX 

lO 

< 



(O w 

in 

O 



Ok «-l 





vU 

^ II 

CM 




•>« X 

K 


« 


o cc o (C 

lO 


1 


O IM X UJ 

(0 


« 


• OK 

• 


1 


K •-! UJ 

CM 


« 


X lAX 

to 


I 


o ^ 

in 


« 


»>4 





Ul 





XO 

o 




tn o cc ^ 

o 




K Ul ui o 

o 




Z • O UJ 

(A 



* 

cc Ul »-• 3 

• 


UJ 

• 

UJ X 1/1 



3 

* 

K UJ 



u 


Ul O cc 




fO 

X 0 3 

o 


X' 


N O VI lA a ^ 

o 


UJ 


X • < »-t o 

o 


> 

J 

Ul O Ul 

m 



uJ 

lA XX !-• 3 

• 



Q 

UJ cc lA 

it 



O 

Ul UiM 




X 

cc ^ X 





O lAK K 

o 


UJ 

ac 

Ul >-• < 

w 



o 

Q K 

+ 

(O 

l~ 

u. 

to VA < UJ 

w 

1 

XI- 


K Z cc -J to 

O' 

> 

O 3 

lA 

CM t-« Ul UJ O O 

K 


Z X 

CC 

• K -J Z <M 

CT 

. tn 

3 ^n 

UJ 

K CC Ul O < 3 

w 


< 

>- 

X X Z <t 

o 

UJ 

-1 Ul 

UJ 

o ^ UJ <z 

CO 

o 

u 

X 

O K Z O 

• 

o 

J <« 

< 

OJ r-# »-tX 

w 

z: 

< 0. 

(C 

• O -J K K 



X irr 

< 

CD o lA <t ^ 3 < 


cc 

cc 

a 

iOOi^i-«tACM X> 


tu 

O X 


o • Ki-«wac*-iuj 


>- 

Z L> 

o 

o u Z _ ^ UJN -J 




UJ 

0> UJ Ul lA K < Ul 1 

o 

J 

UJ 3 

1- 

K lAK < 

o 

»-4 

-1 < 

< 

(AW O -J -13 3 O 

o 

1- 

H- -J 

3 

w ZQ.U14- ZZ CM 

o 

_J 

K 

3 

w O U- X O 

o 


3 J 

O 

Ul (A KO-JOS* K O 

o 

s 

X <5 

J 

CC wzz X UJ_ ^ o 

o 


in X 

c 

'UJ OUJ QXUJUJ z 

• 

o 

cc 

o 

KlAwi-«XOOXX UJ 

o 

LL 

UJ o 


UJ0CK3KXKKK CC 


Ji 

O 2 

» 

X ol «X w ^ 


:e 

< 

I 

KN|t^3PXOU.U. tA 

o 

N. 

OL 

» 

Z UJ W K C0O O 

w 

< 

1/1 


WXJU.UJW UJ 

+ 

U1 



WOIA3BUIZZ O 

to 

< 

Ul 


KZCQ 3 XOO CC 

o 

2 

J 


. Xw<K UJKwt-i 3 

it 


X «> 


O KZXJ) KK O 

K 

2 

X 1 


*-• K Jl Ul O 3 K < <I vA O 


o 

-< # 


UJ X W U. <»-*w O 

CM 

w 

X 1 


XOOO CC >> 1 

it 

►- 

Ul » 


w3<QCOKUJUJ 

• 

3 

1 


OUOCCUIU.X33 

w 

H 

» 


3 X J OK O 


«-< 

«- 


o O <C C w 3 3 


ac 

1 


-13 -j-J UJ r 

O 


» 


030<t KX<< 

o 


1 


O K o X < 3 0 

+ 

►H 

» 


3 -J wo JK Z Z 

■o 

3 

1 


Ul O aJ K K 

K 


» 


K X X K a. <K K 

CO 

U1 



IA3WUJ00 31A01 

K 

o 



3 -«XK> 3 K Z . J 

W 

oc 



3 UJ 3 11 O 

(^ 

3 



3X t 1 1 1 O 1 X 

o 

o 



< 010. 0. 

• 

(A 



IXMKQ- <uI(A 

cO 




1 3300IAOX 


Q 



X<Q-3>>- wwj 


2 



XMKQ3 3MIA(AK 

O 

< 



Ul X — 

o 




X lUO. 

o 

111 



<r M K o 

» 

1/1 



X <K 

CM 

»-• 



< J 

to 

o: 



X 

it 





w 

Q 



X X 


3 



UJ«iJ • 


o 



K K O 

w 




•c < z 


u 



-1 J 



D-10 


D.2 


NASA/MS FC MULTILAYER DIFFUSION MODELS PROGRAM EXAMPLE 
OUTPUT 


This section contains a listing of the output produced by the NASA/MSFC 
Multilayer Diffusion Models Program using the example input data shown in- Figure 
B-1, Appendix B and that portion of the data deck produced by the Preprocessor 
Program for HCl concentrations for a normal launch of a Space Shuttle vehicle. The 
input data shown in Figure B-1 instructed the Multilayer Diffusion Models Program 
to use Models 3 and 4 to calculate maximum centerline HCl concentrations, dosages 
and time-mean concentrations and to calculate the corresponding isopleths. 

As shown below, the Program first lists a title page containing the data set 
title, the meteorological case, date and time, and the height, range and azimuth 
bearing of the stabilized cloud. Next, the Program prints the data parameters for 
each layer and the results of the parameter calculations (for example, mean wind 
speed, wind speed shears, etc.) in each layer. The inputs for Layer 1 are shown as 
the second page of the example listing. This page is followed by a list of maximum 
centerline concentrations, dosages, time-mean concentrations, the time of cloud 
■passage and average concentration at various distances (ranges) along the cloud 
trajectory (azimuth bearing). The Program then produces a series of pages con- 
taining printer plots on log-log scales of the maximum centerline concentrations, 
dosages and time-mean concentrations at selected distances along the cloud trajec- 
tory. 


After the printer plots are generated, the Program lists the range and 
azimuth bearing from the launch pad to selected isopleths of concentration, dosage 
and time-mean concentration. The printed ranges and azimuth bearings are ordered 
such that the output values begin at a range close to the launch pad and move in an 
azimuth direction clockwise around the isopleth. The listing for the Model 4 cal- 
culations ends with a printout of model input parameters for the layers that con- 
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tributed to the ground-level concentrations and dosages previously listed. The com- 
puter Program then lists similar output for the Model 3 calculations except that 
model input parameters are given only for the sample layer used in the calculations. 

The NASA/MSFC Multilayer Diffusion Models Program also contains pro- 
vision for generating plots of maximum centerline concentrations, dosages, time- 
mean concentrations and isopleths of these parameters. Figures D-1 and D-2/ respec- 
tively show computer plots of the Model 3 and Model 4 calculations of maximum center - 
line concentrations for the example case. Figures D-3 and D-4 respectively show 
plots of Model 3 and Model 4 HCl dosage isopleths for the example case. 
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FIGURE D-2. 


Maximum centerline HCL concentration from Model 4 generated 
by the plot routines for the NASA/MSFC Multilayer Diffusion 
Models Program. 
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FIGURE D-3, Isopleths of HCL dosage for Model 3 generated by the plot routines 
for the NASA/MSFC Multilayer Diffusion Models Program. 
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FIGURE D-4. Isopleths of HCL dosage for Model 4 generaged by the plot routines 
for the NASA/MSFC Multilayer Diffusion Models Program. 



APPENDIX E 


DERIVATION OF THE VERTICAL TERM APPEARING 
IN EQUATION (3-18), SECTION 3 

The Vertical Term defined by the summation in Equation (3-18), Section 3 

of the main body of this report was derived using the method of image sources 

(Slade, 1968, p. 346) to account for the reflection of gases and aerosols at the 

earth* s surface and at the bases of elevated inversions. Consider a point source 

located at a height H above the surface (z = 0) within a surface mixing layer of depth 

H^, as shown in Figure E-1. The numbers appearing at the hei|fht z in Figure £-1 

denote the intersections of rays from the real and image sources with the height z 

at various distances downwind from the source. Thus, the number 1 indicates the 

ray intersection from the real source at height H with the level z and the munber 4 

indicates the intersection of the ray from the image source at height (3H^ M - z). 

The heights of other image sources above the z plane are indicated at the left of 

Figure E-1. The parameter yj.is the fraction of material reflected at the earth’ s 

surface. For complete reflection, <^is set equal to unity while, for no reflecticm, 

is zero. Note that the height term for each image source is multiplied by '^for each 

ray intersection with the earth* s suirface that occurs prior to the intersection of the 

ray with height z. Therefore, the image source at the height (2H^+ H - z) is 

multiplied by ‘^^for the reflection at the earth* s surface, from the image source at 

2 

(H + z); the image source at (2H^ + H + z) is multiplied by for the reflection at tiie 
earth* s surface from the image sources at (H + z) and (2H^j^ + H + z). The contribu- 
tions to the vertical term in Equation (3-18) from the real and image sources, 
assuming a Gaussian distribution of material in the vertical, are seen from an 
inspection of Figure E-1 to be given by the following expressions: 
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FIGURE E-1. Schematic diagram of image source configurations. 
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If we sum over all sources and combine terms, the total Vertical Term becomes 

21 f / vSTl’ 
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Where, for convenience in writing Equation (E-1), the quantity 0° is defined to be 
unity. 

Next, consider a vertical line source of length L centered about height H in 
Figure E-1. Assuming that the line source is comprised of an infinite number of 
point sources, the Vertical Tferm for the real line source V^^ centered at height H is 
given by the ejq)ression 
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If the substitution 


? = 


z* - z 

•/i 


is made in Equation (E-2), the resxilting expression is 
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° A / 
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Using the definition of error functions (Abramowitz and Stegun, 1964* p. 297), 
Equation (E-3) can be written in the form 
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Inspection of Euation (E-4) shows that the numerators of the error function argu- 
ments e:q>ress, respectively, the distance of the top and base of the line source 
from the height z. The corresponding vertical terms for point and line sources in 
the surface mixing layer, including provision for partial reflection at the earth's 
surface, are given in Table E-1. 


The Vertical Term of Equation (3-18) can be obtained from the line source 
terms in Table E-1 through use of the following relationships: 
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TABLE E-1 

CORRESPONDING VERTICAL TERMS OF A POINT SOURCE AND LINE 
SOURCE IN THE SURFACE MIXING LAYER 
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where Equations (E-5) through (E- 8 ) ejq>ress the relationship between the height 
coordinate system of Figure E-1 and the generalized height coordinate system used 
in developing Equation (3 -IS), 

Substituting Equations (E-S) through (E- 8 ) into the line source terms in the 
right-hand column of Table E-1 and collecting terms results in the expression 
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where, agaia 0** is defined equal to unity. In Equation (3-18), the factor 1/2 
appearing before the bracket in Equation (E-9) is contained in the form factor 1/2 
■^2 it appealring in the denominator of the first term on the right following the equal 


